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I 16 Ab^pect 

This computer user’s manual describes the operation and the essential features of the Microphone Location 
Program, one of the two programs developed under the Method of Fan Sound Mode Structure Determina- 
tion Program, NAS3-20047. Jointly, the two programs are used to determine the coherent modal structure 
in annular geometries. The purpose of the Microphone Location Program is to determine microphone 
locations that ensure accurate and stable results from the equation system used to calculate modal struc- 
tures in the second of the two programs. As part of the computational procedure for the Microphone 
Location Program, a first-order measure of the stability of the equation system is indicated by a matrix 
"conditioning” number. 
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Ml mni) Ol I AN SOUND MODI SI lUK I OKI Dl 1 1 UMINAI ION 
( OMI’llI I U I'KOOUAM USl U‘S MANUM 
M!( KOniONI IO( AIION l’K(U;UAM 

In 

(i. I . IVkiMt. K. A UHls ;mtl k. A I ,m> 

\ Whitiu’v Aiivi;itl Omiip 

l.n SUMMAkV 

This lominiU'i um’i s mami;il ilosi iilvs llu' i>iH‘ralioii .itui llio osM>tiii:il IV.iIiiivn of tiu' Mivio 
phoiii' I oiMlion I’lOiiiaii'. otic of llic two I'ro.i’.r.mis i!cvcli>pcii utuici the Mclhoii of I an 
Soiiiul Miulc Slnictiiiv Dcicriiiinalion I’rojiiam, NASo-.'lU)’} .KmiiiIx. the two pi'oiu'anis ate 
usci! lo liclcfiniiu* till' coliciviil iiu'dat slmcliiiv in aniuilar lici'inclncs. I lic piirpiisc iil the 
Microplionc t ocation i’ropaam is lo dcicnninc microplunic kicatioiis ihal ciisiiic acciiralc 
aiul slahlc rcsnlis lioin llic ci|iialion s\slcm used lo calcnlalc iiunial slmcliiivs in ihc second 
ol llic iwo piocaams. As pail of llic coinpnialional procedme lot llic Microphone I iication 
I’focaam, a lirsl-otdcr nicasuic ol llic stabiliU of llic ciinaliim sxsicni is iiulii.alcd In a nial- 
ri\ "condilioiiiiii*.” iinmbcr. 


:.0 IN I kODl 'CHON 

New fan dcsic.tis for niodcrn liich b\ pass ratio commercial cncincs iilili/c biadc-xaiic inicr 
aclioii ihcotA to the cMctit possible for contiollin!’ the pn'pae.alion of inieiaclion noise I'm- 
reniK, this llieory detines the modes fh.il can propaeale, but has not been rieveioped lo the 
v'Xtciil Ihal il can leiiably predict Ihe slrene.llis ol the pia^paealin;', modes. 

binlher luiise rediiclioii vamiKI be achiered if Ihc piopat’.aliiie. modal stnicliire were qiiaiili 
lied. Once Ihe imulal slrnctmc were defmed. an aiiab lical '^vslem t\>r acoustic Irealmeni de- 
sii’.ii could be iilili/ed lo opiimi/e trealmeiil for a jiireii modal sliucinre. to produce more 
efficieiil schemes In addition. Ihe modal structure could be employed lr»\eiif\ ileveiopiii.e. 
theories of fan noise iieneralion. l o provide this cai\ibilil\ b> means vd' measured data. Ihe 
Method rd I an Souiul Moile .Siruetnie Delerminalion I’ronram (N.A.S.t .’DO-I I w.is imdert.iken, 
1 111’ melhird would be ulili/e«l iinld ,i x.ihd l.m noise I'.encr.ilion nu>dcl on .1 iiioilal basis be- 
comes a\ ailable. 

Ihe Iheorx upim which l.m spmniiii'. im>dc Iheoix is tonnded w.is piesenled m D’nl b\ I \ ler 
•ind Sofrin del 1 I. followm.!’ exieiisne aii.iK li mI .iiuI evi'enment.il sliniies 1 ,ilcr. Sofim .nui 
Mi( .11111 del .’t derixed Ihe I'.eiieial lorm of .1 crdieieni acoiisiiv- xx.ixe m an infinilelx lone, ex 
hndiic.ihiuci w Inch cNtended Ihc Iheoix 10 mchide elfccis ol .i\ial flow I his e^iualion ex 
ptesses Ihe cohereiil .icoiislie presstne ,it loe.iliims in the duel .is .1 Inm lioii id ihe .nn)dilude 
.md phase td Ihe pi op.p'.iiin;’, modes eompiisine, the sihiiuI lield I hese puieix crdieienl sii'. 
n.ils, w hteh .lie due lii Ihe eonltibnlions id ihe consiiiiieni modes, .ne exti.icled fiom Ihe 
oxer.tll sien.il in enhaneemenl leehniiiues .id.ipted .il l*i.ilt \ Whilnex Anci.itl Ihe .idx.m 
l.iees ol utih/ine sij’iul enh.uicemeill is discussed In I’lnex m leleienee ' 
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.< 0 l’KO(;UAM Dl sSi KirnON 


.1.1 AUiOKIIMM 


— 

I ) I lu‘ iminlvr ol M li'k iod KkmMoiis i-i|imls lIu' mmilvr of tiunk-s, 

svl^lur'"' M .„o,l „ „„ a„ 

kv, .i„- 

- riK> kv..lions .,V .vsnik tk-J lo tiu- ,Uk., w,li u,i(il liu- sivalu' nunl,s lv...nu. 
Pi.ulK.ill> iiulishiuimshsihL- i!u>n p,oK- KuMlions jr. usoJ 

‘ or Iho SRUIUI !u-kl ,n (Ik- JuvI .omni .M,,,. v 

tliK t inoiK-s, (Ik- i'.i-oiiiv'Itik- i>.ii.iiiK-(,.|s ,, i,, .t , , i. ■ i I ) N.kounIk- 
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N nuuik's .uv llu'ii i»nli‘ivil siu li lli.il Ilk' k iivunilViviili.ii iiknU' iikik'os .iiv ;iir.iii|k'il l.iijvsi |i> 
Mii.ilU-Ni .Hut ilu' k-i'iu"<iumiiiiit’ r.iili.il iiulii-ios .nv nuU'ivil Horn laiik'Nl lo sm.illi'sl I'luli 

iiii-iimlViviUi.ll ilUkti' niik-v 

rik’ iiiK iopluuk- lik'.ilioii piik'i-ihm- i^ iuili.ilikl b\ inliiklik iii!’. I'otli Ilu- liisl iiuuii- in llu- 
stk’> ilK-ii oiili'i .iiul Ilk- alkUl’ ivtViviia- iiiii-rupluiik- ku-alioii iiilo tlu- vi|iialioii s\si(,-m 

I lk- M-n>iul luoik- 111 Ilk- si-iiiK-ntial onk-r is ilu-n inluHlui-i-il inio ilu- t-i|ii.itioii s\sli-ni aiul llu- 
1.-01 i\-spoiuliii!’. si-i-oikl Ku-alion is ik-U'iiiiiiUki using a stoi-lu'slu- sv.iicli puk-i-iiiiri-. I Ins piu 
(.-i-Jiiiv sk-k-i-is a kk-aliim llial ma\iini/os Ilk- absoliiU- magnKiuk- ol' Ilu- ik-U-rmiiiaiil assiuiaivtl 
wilh Ilk' ik'siilling ' \ maln\ oipiation svsU-in. Tlu* si'li-fk-i! iik-alu>ii is oblaiiu-il b> an in-ia- 
ii\i- puvv-ss in Ilk- iko i-ikn\linalos on tlu.’ iliu'l wall using laiulonib gi-iu-iau-il kualions 
in Ilk- ivgion klu-iv niivrophoiu's can lu' plaa'il. 

rtk' Ins! ailt-nipl gi-iu-ralos inic ioplu»nc- cooiilinaU-s t'l oin a tlaussian ilisirilniiion c'liarac ii-i i/i-ci 
b\ .1 nu-an al tlu- niiildK’ ol ilu- sc-aivli ivgion wilh a staiulaul ik-\ialion c'c|iial lo hall' Ilk- k n- 
glh of ilk- soaivli ivgion I'ivc- tunulivil caiulidali' tnicioplioiu- lovaiions aiv gi’iu’iau-il wiiliin 
llns ivgum. aiki tlu- c-onvspoiuling ik-U-nninanl of tlu- c-ipialion sssii’in i.s caUulaU-il. Only llu- 
lhin\ micii'i'lioik' lovaiions lhal \ ivk! ihv largi-st inagniluck- ol lliv cU-lvnuinanl and Ihv vaiiiv 
of Ilk- dvU-miinanl aiv iviainvd. Tlu- nu-an and slaiulaui dc-vialuiii ol llu-sv ihiilv valiivs aiv 
valcidalc-d and llk-n usc-d lo ivsirivl ilu- lU'M altvmpl lo dvlvuniiu’ a suilablv lovalion. 

I 'lw hiiiulivd micTciphoiii' tov.Uions aiv lunc raiuloinh gviu'ialvd lor a sc-vond tiinv. bill voi- 
ic-spondmg lo Ihv lu-w stalislival inUn nulion. I’liis piocvss of ivdc-rining Ilu* stalislivai para- 
nu-ic'is .iiid I'.indoinb gvnviating tiw luindivd lovalk>ns is vcniliiuii-d iiniil itiv langv ol ihv 
lliiil\ c.dik-s voi lv^pondlng >o vavli mivioplionv voc'i'diiialv has vonwigvd !■.> wilhin a spvvtik-d 
Uik’iaiKv or uniil tifiy ilvrations haw bvvn voinplvU'd 

Al ihis link'. Ilk- svlvvlvd iniviophonv lovalion wliivh vonvsponds lo ilu- iaigvsl vahiv of 
ilk- dvivrmiiiant lalviilalc-d licnn ihv ivsulliiig ’ \ - ni.itii'c viiiialion syslvni is v\aiiiiiu-d 
In Ihv niiviophonv lovaliiui avci-plaik-v piovvdinv. voiulilioiiing miiiibvi iiv!'. >) is c.ilvn- 
lalvd as Ilk- lalio ol' thv i-\liviik- vigvn valuvs assovialvd wiih ilu- vi(iialioii syslvni. I lk- von 
dilioniiig iiiinilvi is a tiisl-oidvi mvasiiiv ol Ihv sviisuivil\ of ilu- vc)iialion syslvni lo small 
pvi liiibalii*ns in mivrophonv lovalion c'r pivssiirv nivasiiivnivnl ll Ihv voiulilioiiing mmibvi 
IS k-ss Ilian Ol vipial lo a spvvilivd vahiv |v. g . oiiv lunulivd), whivh is inpiil lo Ihv piogiain. 
tlk-ii ihv vipi.ilion s\slvm is vonsidvrvd U' bv wvll i-ondilimu-d in .1 in.ilhvinalival svnsi-. In 
llns iMsv. Ihv inivioi’lionv kualii'ii is avvvpiablv, and Ihv llmd mode- in ihv ordvivd svc|uvik v 
is iiiHculiivvd III Ihv ivsiilliiig, } \ vi|ii.ilion s> stvm. 1 liv sli’ichioinvltiv svaivh piovvdmv is 
llk-n miliaUd lo kualv a ihiid inivuiphoiiv mi Ihvdiivl wall, using a two I'aiamvlvt tX.fU 
svaivh pi ovvdin V. 

Whvii Ihv vondilioiung iminlvt is giv.ilvi than ihv spvvirivd v.ihiv, luo tkissibililivs vsisi In 
Ilk- liisi possibiliiv. Ihv sioilu‘sitv svaivh piovi-dmv is iviiislalvd k' c>blain a suilablv nuvio 
plk'ti. k>v.iiii>n on a udi.il pu*bv l lns proivdmv iisvs a two paiamvlvt ii. (I) svaivh, assiiimng 
,1 i.idi.d pu'lv at ihv rmllk-si lov.itu'n liom Ihv fan kilhin Ihv ,i\i.il ivi'U'ii Siibsvciiivnl niivio 
plionvs ,iiv pl.iv, d on ilm piobv iisiiig, .1 onv-p.ii ,imvk-i ii I sv.iu h mill! ilu- c oiulitioimig mini 
bvi IS |-.iv,ik 1 lhan ihv spvvilivd \ahiv 01 ihv m.iMiiiiini miinbvi ol iiik iophonvs pvi piobv is 


\ 


•suttsliai. In niluT iwn iMnnuMoi (i. 0 ) s.iiivh is use,! In seKvl ;in altem;.le loenli.n. 

. ssuimni- ,t nev' i:uli;)l pinho. i |u- smnul pnssihilily is wlien (he sek-t ie,! miernplioiu' l,u;i 
lloiMs IIK' lust in Iv i.liiaul nri pmho, using tlie iibove lwo-par;imoler (r. 0 ) soardi aiul tiu« 
comliliomng iiumK-r is still givatn llian Hie spccirial value. In lliis ease, the Ineatinn is ae- 
eepled aiitl no lurllier iniemplmnes are permitlcl on the piolv. Kadial probes are used to 
oeale miemplmnes loi eillier possibility until the maximum number of probes as speeilied 

>V input IS leaehed. When this oeeurs, a message is outpiil to the ellVel and the program exe- 
eulion IS termnuled. 

As an ,.ption, the Mieroplione I oeation I’rogram also ean be used to a.ssess the suitability of 
an eMstingsel ol mieroplione loealions to provide salisfaetory input data to the Modal (‘al- 
uilalion I log, am, I lie previously discussed stochoslic search procedure and acccidance pro- 
cediiie are not uli ,/cd m this option. The .suitability is a.ssessed by means of the value of the 
conditioning number calculated Irom the eigen values of the ecpialiou system. If the number 

is less than a specilied level, the c,|ualion system is considered to be well conditioned in a 
mathematical sense. 

lixperience pained from running analytical test cases has shown that a good uptier bound for 
the amd.tionmg number is on the order of ten. Iu>r this conditioning number, at least three 
MgmtKjnt tiguies are retained in tiie matrix inversion algorithm. 'I'o expedite the intermed- 
iate steps ol locating each successive microphone, the value that is internally initialized at 
execution ol the program in lieu of any input value has been cho.sen to be 100. This value 
can be altered by input prior to the data case if more or less accuracy is desired. I lowevei . 
to uisuie stability, the linal conditioning number should be about ten. 

3.2 1'R0(;RAM OVl-RVIL’W 

Ihe Microphone I oeation Program compri.ses six program .sections which are utilized in part 

or whole to accomplish the objectives of the two possible modes of operation. These six sec- 
tions are: ■ I ^ 


1 ) Input - The input of all data is by the NAMlil 1ST specification, ami tlie internal para- 
meteis aie initialized lor program execution. 

2) Characteristic Number Calculation - The characteristic numbers k " and O" are 

calculated using Ihe procedure in Appendix A. ^ 

.^) Mieropimiie Location Determination An iterative procedure is defined tor siiccessivelv 
locating each microphone. I'tie algorithm tlelermiiies loealions that maximize Ihe de- 
terminant Ol an e,p,alion .system characterized bs perturbing one location while Ihe 
piivious loealions are li.xed. I his procedure was tiescribetl in Section ,L I. 

■D ( onditioning Number Calculation - A conditioning iiimiber Iref, is calculated fm the 
acoustic wave equation matrix after delerminatioi, of each microphone location The 
conditioning number is .lefined to be the r.Mio of the extreme eigen values ,.| Ihe matrix. 

M Microphone i oeation Acceptance AccepMnee ,»f a micophonc loc.iliou is b.ised o„ Ihe 
iiiagmiiideol ihccondili miiigimmlvr. Ihe location imucdiiie is uiiisl.iied v.hcn the 
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viimlitionim- mimivi is )>,iv;iU'f Ih.m ;i sTil^ mini siH-i'ilivd ,is in|uil. A v;ilm' o\ HU) is iii- 
lonnilU il ii \ahK‘ is luM siiin'lii->l I'V U^' A lI.-su ii>tioii ol tlu' minoplioiu' 

;ut\'plaiKv proivtimv was pivsnili-ii in Soilion .VI. 


Ml Oulpiil - All ivsulls mlailat^'il by tlu’ pioinam aiv piinlinl. 

Hu- inii'nvlationsliips Ivlxuvn llu- sis pmi'.iam M\lions ami llk'ii ntilils lor isu-li opium is 
illiisliaUki in l ij'.uiv 1. As iiipni. luitli options ivipiiiv a sp.k-iliv' imuiv uronik ink't jn-onuMry. 
ami U'sl anuiilion to valoiitaU' cliaracK'iislic mnnlvts. Om- oinion. A. ivr|nnvs midilional 
iiunii in tlu' lonn of suvliastiL- s.-aivli iniranuMm-s lo lummi llm I'.mmicliA wlu-iv nmnop lom- 
kH-alioiis oan lx- i>laa‘il. Aililitionallv. oplion A osainim's tlu' anuiitionin}- mimlu'i to iliMoi- 
mi,u' w lK-lluM- IlK' minativMi syslmn assodatorl willi iho introiinction of mull snavssnv inodc- 
mimopliono lomtion is svoll .omiilionmi in a inatlu'nialiml sonsi-. l’lu> otlu-r option IV ivqimvs 
tlial tlio inicrophoiu' locations Iv spccilVd as inpnl citlK-r arbitranly or as a scl ol csislinit 
locatimis In both oi'lions. tlio comlitionin}’ mnnbor is calculatml lor liio acoustic wave crpia- 
tion syslcin. Tliis mnnbcr is nlili/cd by tile coniimtcr inoitram in option A to dclcmnno wlictlicr 
■■acli microphone loealion is acceptable lor acriniiinp. esivrimenlal .lata. Option H supplies 
this value to the user to evaluate an esislinp. set of locations. I’lie results Ironi both options 
are i>rinteil b\ the output section. 


.V.^ IMUXIRAM SUHUOU riNl S AND I'UNC TIONS DIVSCKII’TION 

riie subroutines ami funclions use.i in ihe sis proitram sections presented in Section .V: are 
listed Ivlow. riie purpose of each subroutine .n- fiuietion is dcNcribed. Also, as approin iale, 
a iM-iiiciple element diagram of the more complicated sections are presented and discussed. 


Input Section 

I lK- input of data to the comiuiter pio.irain is In the NAMI I IS I format. I his form of in- 
pul is described in Section V4. 1 . and are set by sivcifyinp botli Ihe mpul variable name and 
its value. In .Section .V4.1 Ihere is a lisline of Ihe input variable names with a correspondini: 
descriplion of their iniri>ose for oivration of the computer pmitiam. All in|ml are read into 
the proi’.rani In Ihe folknvine, snbrouliiie. 

INI’Vr - riiis subroutine in|itils data for each case and sets up Ihe ueccssaiv inleinai paia- 
meteiN. 


(.'haracteristie Nuinher ralctilation Section 


i Npiessions are derived in Appendix A for solviiiji two simullaneoiis ciiiations that vield Ihe 
diaracleristic numbers k; " and 0'’,^. A principle elemeni vlianram is ,nov ule. in l■.!■.ule ^ 
to ilhMiale Ihe funclion.il elements that lead lo a delermmalion ol these numbers. Imtiallv . 
the Older of Ihe Hes^el limclion ,1m and ^ ni are determined from the circmnlereiilial ordei 
of parlicular mode. 1 he .lin and Vin Iknvl funclions are then evaluated, as appropriate vie 
pcmlms’ vm Ihe v.thie of the duel hid' tip lalio v ieldiiu' Ihe value ot and 0,,,^, bv solvnie 
tlie simull.meous eviuatii'iis eomptiMtie these Hmsel limctiv'in 
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1 ^. 


















1 III’ Mihrniitiiii'^ iind riinclioiis l|);ii ;iiv iilili/rd in lliis si-iiinti ;iio 

Kl.H'AI - Tliis siiliioiilim' inkiil;iti‘s llio di;ii;icteiistic numhi'is k!", inn) ( > " 

HI/' • inp 

KMUCAI - Tlii?: siikrouliiu' is iisi-il hy K(.K'AI. to i-aliu!;iU‘ llio i‘hai:u.Mi>iislic mimlior 

I MMCAI - This Mihrmiiitu' lak'uliili's i hiirai-lnislic 1 iiitu-lion vaUu-s |oi :i piirtii-nlai nulial 
valuo r/lv 


l•■Al./.ll’ - I hi-.. riiiiLtit)n soivi's loi a root of a pivon riiiuMion iisinfj a L'njiihinalion of false posi- 
tion ami hisi.’ction iedini(|iies. 

HHSI.l - This function is useii hy KMUC'AL to calculate values of for the ci|ualion 
wiiich itefincs the .system of dilTcrcntial equations. 


dr ' 


lir 'L 
d 




- Mn,(ak- )|+Q« ~ 
dv' dr' 




= 0 


= 0 


for a luih-lip ratio not equal to /.ero. 


BRSL2 - This function is used by KMLK’AL to calculate values for the ecpiatiun which 
defines the above system of differential equations for a hub-tip ratio equal to /eio. 

BHSJ - This subroutine calculates values of the Bessel function of the first kind. 


BUSY - This subroutine calculates values of the Bessel function of the second kind. 
Microphone Location Deterinination Section 

The algorithm for ileterniining microphone locations was discussed in Section I and illus- 
trated by a principle-element diagram in I'igure .V Initially, statistical information is calcula- 
ted for the geometric search range l>y determining the mean and standard deviations of the 
axial, circumferential, and radial axis where microphones can be placed. Llie next succcs.sive 
mi''rophonc localk)ii is detei mined by restricting alt previously determined locations. I his 
location is obtained by generating tlve hundred locations and retaining only the values that 
yield the thirty largest determinates of the erpiation system. If the tolerance of the thirty 
values associated with each microphoi,c coordinate is witlim the specific value, a .suitable 
microphone location is determincil. Wnen the tolerance is greater than the specilied value, 
the search is reinstated tor at least fifty times by adjusting (he statistical intormation to re 
strict (lie next search in the geometry associated with the hugest determinates. The sub- 
routines and functions that arc ulili/ed in this algorithm arc: 

S roCIl - i his sidnoiitinc performs a stochastic search to determine microphone localions. 
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INIilA! I iiis siihnuiliiu' i-. liv'd hy SrOCll lo initiali/v liu- lU-ivssary parainoli'is loi llii- 
slni'haslif si'aivli. 

S TM ■ Tliis Mihroiitiiu.’ is iisai hy Si nCII to iviinaU' tiu' various statistics iv(|iiircii for the 
stocliaslic scarcli. 

random - I his suhroulinc is used hy S I’OCH lo {icucialc llic random values associated with 
tlie microphone locations useii in an evaluation. 

S rVAR - This sidnouline is used hy RANDOM lo provide random variables lor either a nor- 
mal or an c.xponenlial distribution. 


RAND - This subroutine is used by STVAR to calculate random numbers. 

DITC'AL- This t'unclion is used by S TOCH to calculate the determinant of the matrix equa- 
tion system. 

UPDA'I'I'. - d’his subroutine is used by SiOC’ll to upilate the table ol ileterminants with the 
current value of the matrix determinant. 

I’RIN'IT - This subroutine is used by STOC'H to print the results of the stochastic search. 

CONVRO - This subroutine is used by STOCIl lo test for the convergence of the microphone 
coordinates. 

Conditioning Number Calculation Section 

The conditioning number (rcl. ,^) is determined as the lunctional elements illustiatc in I ig- 
lire 4. Initially, (he eigen values ol the equation system tire calculated and ordered numeric- 
ally. In this way, the ratio of the maximum and minimum eigen values can be computed. 
The subroutines that calculate the eigen values are listed below, 

lilCCC, I'HALAC, i:ni:s.sc. IXRlllC, l-l.Rn’t'. 1;B13CKC. U1:RTI SI - These subroutines 
calculate the eigen values and eigen vectors of the equation .system from which the condi- 
tion number is calculated. I his package is the property of International Mathematical and 
Statistical libraries, Inc. of Houston. Vexas, 


Microphone l.ocafiou Acceptance Section 

In Section .V 1 was a discussion of the procedure for determining whether the selected micro- 
phone locations insure a numerically stable solution when used in the MF( . A principal cle- 
ment diagram is providi'd in I'igure 5 to illustrate the acceptance procedure. I he logic en- 
closed by dashed lines represents the acceptance procedure. I he inicrrelationships between 
Ihis procedure and other program seclioiis is illustraleil by the log.ic lines extending outside 
of the dashed lines, luitiallv. the conditioning niimber is compared lo the acceptance value. 

If the conditioning number is acceptable, Ihe microphone localion |irocedure is reinitiated lot 
the next mode-microphone combination. If the conditioning numbet is (ucatcr than Ihe 
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;u’a’|i|.iiui’ Viiliii . Ilu' miv roplioiu' loialinii i>ioivdiiiv is ivinsl:Ui'il In m-KvI ;in ;ilU’in.i;i‘ In 
i Mlmn on a [iiolv. Ailililioiialh , (lu- alUnvabU* liinil loi mii rophoiu's |vi piolv ami mimlvi 
ol piolvs au‘ i lu'okoil Sina* this pioivihiiv is |vi romu’ii in tin; MAIN, ilu'iv aiv no snlirou 
linos or Inrniions li> list. 

Ontpni Soo(i*>n 

I III- oulpiil lormal ami variables iVoin oporalinp, llu* Ml I’ aiv ilisinssal in Soolion .V5.1. A 
sample ease is pre.senleil in Seelion .V.S.2 lt> ilhisirale the exeeiition of a ease eompi'isinjt 
three proinitiatinfr nuHles. Both seetions mlrhvss the two possible imnles of operation tlial are 
e\eeulahle with this proj'.ram. Results from the eomputalional proeeihire are printed by the 
subroutine lister! below after all antiles are ei>iiverteil to within the range 0" to dot)''. 

I’RIN I I'his subroutine prints input and output values. 

ANtil’O.S I Ids subroutine eonverts negative angles to positive angles in the range 0" to 
ooO" for printing. 

.1.4 INIHU Db.SCRimON 

.1.4.1 input rormal 

The N AMI I l.Sf format is useit tr> input riala into the Mierophone I oealion Program and 
eonsists of a list r>f parameter names grouivii muler an iiientifying name: AINHA I’.A. The 
parameter names eonespruid It' variables single variables ami matriv elements used in 
the proi'.ram. These variables are set by speeilying both the parameter name and its value. 

■A feature of this type of input is that all a.ssoeiated parameters need not be speeifierl. Any 
parameter not speeified in the i;ipul retains its value irom the preeeeding ease or the default 
\alue if the input is for the first ease. 

N.\MT 1 I.S T injnil for eaeh ease is iilenlified by the eharaeleris &INn.A T.A in (.‘oliimns J-7 
of (he first input eard. Reghtiiing in ('ohnnn parameters may be set using the formal: 

Parameter Name t'onstant 

I he eoiistanl may be either a real or integer value and must be followeil iminedialely by a 
eomma Parameter names, assigned values, or neees.sary eommas must not extend beyomi 
I ohmm ' ami names or values eannol be eo itinued on a subseipient eard. imbedded blanks 
are not ivrmilteil in either the paramelei name or eonstani value. Paiameter names ami 
(hen assoeialed \aliies may be speeilied in any older. The eharaelers tATNI) signify the end 
ol the in|nil lor a part ienlar ease. If .uidilional eards are reipiired, paiameter names must 
begin in ( ohmm \ 

■A s.nnple ol ihis lorm ol input loi loealing, three mieiophones is piesenlerl in Tig.ure <>. 
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.<.4.- Input 1’itranu‘lcrs 


A sign convcntiori was aiioptcd for assigning positive or negative values to the input paranie' 
ters. Any input parameter not aiUlresseii in this cliscussir)n is a positive value. The sign con- 
vention is tormulated witli respect to a cyliiulrical ciKrrilinate system tiiat is consistent with 
the derivation of t!ie coherent acoustic wave propagation model. Its unit vectors are designa- 
ted hy the directions: axial - x, circumferential - 0 , and radial - r. 

A constant radius, annular duct is aligned with respect to this coordinate .system so that the 
positive axial unit vector projects in a direction opposite to the flow. Thus, the Mach number 
of a uniform axial flow is always designated by a negative value, denoting the axial flow rate 
in the negative axial direction. A positive circumferential unit vector projects in the direction 
that the rotor spins with negative circumferential values related to the counterrotating rotor 
direction. I'inally. the radial axis projects perpendicular to the centerline of the duct; thus 
radial values are always positive. 

Hach mode is characterized by three parameters which represent the circumferential and ra- 
dial jiressurc di.stribution and its propagation direction. A specillc mode is uniquely defined 
by the parenthetical notation (M. p). The M defines a periodic circumferential pressure dis- 
tribution with M number of lobes. Positive integers represent a corrotating M-circumferential 
lobe pattern with respect to the rotor direction and negative M integers refer to counterrotat- 
ing modes. The radial mode index /i corresponds to the radial pressure distribution. These 
values are always non-negative integer numbers, with high integer values indicating large pres- 
sure variations witli respect to the radius. 

The modal propagation direction in an inlet or discharge duet can be either an incident wave 
propagating from the fan ora renected wave propagating towards the fan. Wave propagation 
in a moving medium is similarly affected by the How rate for modes that are propagating 
wiiii or against the flow direction. Hence, the input variable lOlR designates wave propaga- 
tion with respect to the flow direction. Positive values denote waves propagating in the op- 
posite direction with respect to the flow, such as incident waves in the inlet duct and re- 
flected waves in the discharge duet. Modes that propagate in the same direction as the flow 
are designated by a negative value for the input parameter IDIR. 

fhe assigning of values to the input parameters will now be considered. 

Since a determinative equation system is required, the number of microphone locations to 
be calculated is equal to the number of modes indicated by the user. To initiate the selection 
proceilure, the input must define at least one microphone location at a reference point. In 
order to reduce the experimental commitment by avoiding probe-mounted microphones, 
the referenced location should be on the duct wall. When option H is utilized, a number of 
fixed microphone coordinates equal to the number > pecified by the input parameter Nl) 
are required. 
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A ik-.sai|>linii of iho input vniiahk-s is proviilcd in l ahic's I. II, ami III ; Table I Ceneial 
Input Parameters; I able II lest (leometry ami ('ondition Input Parameleis; Table III 
.Sioeliaslic Search Input Parameters. Under the column heading! "Variable Type”; the letter 
"U ’ imiicates that the number is real ami contains a decimal point; the letter "I” indicates 
the number is an inteiicr and docs not have a ilccimal point. “Delaull Values” are also deli- 
neated and indicate the value of the parameter that is internally initialized prior to the pro- 
tuam execution. Parameters not specified in the input for the first case retain this value. 
Although the detaull values are ex|>resseii in units of tlie Ibiglish System, the computer pro- 
giam can be executcil with ilata in any consistent system of units. 
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page K 

^ ntoit QVAurv 


(il-Nl RAI INPUTI'AKAMI n KS 

Input 

Variable 

Default 


Name 

Type 

Value 

Description 

NI.OC 

1 

> 

Niiinhei ol micruphimcs or modes. (Less 
Ilian or erpial to filly). 

II. MU 

1 

[) 

t’riiil iiidieator lor eliaractcrislic L-fuiietion 
value. 

0 = No Print 

1 - Print 

M(l) 

1 


Circumrercnlial mode index. 

MC) 

M(.U 


• 2 
0 

(Inpul NLOC values) 


,M| 50) 

0 


MUS(I) 1 

0 

Radial mode index. 

MLISC) 

1 

(Inpm NLOC values) 

MUS(.?) 

0 



MUS(50) 

IDIR(I) 1 

0 

1 

Mode propagation direetion indicator. 

IDIRC) 

1 

(Input NLOC values) 

11)1 R(.t) 

0 

1 - opposite flow direclion 



1 = with (low direction 

IDIR(50) 

0 



17 


I 


tahi I- II 1 1 sruioM! Tin’ ANnroNPiiioN iNrnri’AKAMi' n us 


Input 

Vaiialilo 

nntaiilt '•'* 


Namn 

Typn 

Viilni! 

Oiisciipiion 

NIK 

R 

0.1 .» 8 

lltih lip ratio 

OR 

R 

.s.o 

Oiiti’i lailiiis ol iliii 1 

1 M\ 

R 

0.07 

.-\\ial Nhuli miiiilH'i (always nosiiivi’i 

IRQ 

R 

.MOO. 

lost (r 0 ;lll 0 lli'\ (Ill'll/’) 

Sl'l i-n 

R 

l.i>t'(i. 

Spooit of souiui 

NO 

1 

1 

Numhor ol iiiiiialh lixoil iiiioiophouc 




kvatiotis (within laiigo ono to Nl IX") 

xNun 

R 

0.0 

■\xial ooouliiiaios of llio I'isoii miciophono 

\M(1) 


0.0 

luoa lions. 

XM(-P 


0.0 

(hipiil NQ i.iliu'sl 

XtvU.SO) 


0.0 


RNUn 

R 

s.o 

Railial oooulinaios of itio fixoil niioiophoti' 

RNU:i 


0.0 

looations. 

RM(a) 


0.0 

(Input NO vahios) 

RNU>m 


0.0 


niM(i) 

R 

0.0 

(■iuiiinl'i'iontial ooouiinali's of iho liM'il 

TiiMi:) 


0.0 

luiotoplioiK- looations. (ilopioos) 

TUMl.n 


0.0 

(Input NQ v.ihiosi 



Oti 


Null*: (a) 

IVlaiili val’.ii's shown in lahh' 

ail' in units ol ilu‘ t'ii};hsh .Snmi’iu. 1 In' ptoj'.i.nii. howi’vi'i. is 


ili'sijMii'il 1i> t'viiiii- with lima ill aiiv omMstonl s\sii'iii ol itml-. 
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V" ■ 1 
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I Aim 111 s roniASTir si aucii ini*ii t i'akami-ti ks 


Input 

Name 

VtiriahIc 

lypo 

i.il 

Default 

Valuo 

Description 

ll’l 

1 

(1 

Prim itulieaint lor laliie of ckMemiimmls 
and mieroplione enoidiiiuies. 




0 “ No prim 

1 - Print |•■ilUll Seardi 

2 = Prim every Seardi 

NSKH 

1 

.St) 

Maximum immber of .seardie.s. 

NVAl. 

1 

soo 

Number ol' evaluations per search. 

IXPNT 

1 

. .2 

lixponent of convergence tolerance, 
(e.g., Tolerance = lomPNT) 

rONDNO 

R 

100. 

Value of the conditioning number for 
inicroiibone location acceptance. 

NPROB 

1 

1 

Maxiimim immber of probes in duct. 

MPROH 

1 

1 

Maxinnmi niimbei of microphones per 
probe. 

XMIN 

R 

0.0 

Minimmn axial location foi microphone 
placeineni. 

XMAX 

R 

110 

Maximum axial location for microphone 
placement. 

TMIN 

R 

0.0 

Minimum angular location for microphone 
]rlacement (degrees). 

IMAX 

R 

.R>0. 

Maximum angular location for miciophone 
placement (degrees). 

mix 

R 

1.0 

Minimum axial disiaiux' between micro- 
phone locations. 

Dlil.R 

R 

1.0 

Minimum radial distance between micro- 
plume locations. 

m L III 

K 

1 I..S 

Minimum angular di.siancc between micro- 
phone locations, (degrees) 


Ntiic: (;i) Ddaiill values sluwii in lahio ate in iiniis ol tlie I iiplisli Sysiem. I lie pmurant. Iiowevei. is 
ilesiniK’il In execiilc with data in any cniisisieni system nf imiis. 
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.V5 Ollil’UT DI SCKimON 


Oiiipul l-ormal 

I ill' iMiti'iii litun i1k‘ Mii ni|>lu>iu' I oiation l*m);rain is ornani/i’il ituo loui sections; Siocliasiic 
Search Ovii|nit, Input Vaiiaolcs. ('alculatcil Micioplurnc t.ocalions, ami Characteristic l -l unc- 
tion. All tour sections arc inchulerl as output when either option is retiuesletl hy the input. 

The computer printout for a sample case is proviilcil in Appemli.v C. 

I he Stochastic Search Output Section comprUes the value ol ihe input parameters that ilefines 
the iterati\e process. These parameters inehule the numher of iterative evaluations, tolerance 
for convergence, ami the numher of microphone coonlinates used to generate statistical in- 
formation. The remaimler of this portion of tlie output includes calculated values from the 
slocliastic search process. These values are the statistical information that corresponds to 
the searcli raiijie for each microphone coordinate and the value of the detenninant associalcil 
with tlie successive introduction of a mode-microphone comhination in the equation system. 
The final ou put value in this section is the conditioning; numher of the equation system char- 
ac'eii/ed hy the microphone location associated with the lar^iest value of the determinants. 

A print indicator is availahle to the user for restricting the depth of output in this section. 
However, ilic final conditioning numher of the equation system is always included in the 
output. 

Input \'ariahlc .Sectio.i includes the value of various input parameters that were supplied hy 
tlie user. The three indices that dellne the modal structure, the circumferential and radial 
i'rdcr. and the wave direction indicator are listed. Under the l est (leometry and (.‘onditions 
Input iieading. there are various parameters that define the fan iluct geometry and operating 
conditions observed during the experimental program. I hese parameters include the duct 
radius, duct liuh-tip ratio, axial Mach miiiiher, and frequency. .Additional stochastic search 
values include the axial, radial, and circumferential search ranges that hound the locations 
wliere microphones can he placed. 

Hie Calculated Microphone 1 ocation Section includes the cylindrical coordinates of the 
selected microphone locations. One microphone will he located for each mode supplied hy 
the user as input. 

Hie final section ( haraeteristic I -iTinctions. includes the value of the I -functions I tk,jjjui t. 
at the measurement and prediction locations corresponding to each mode in the sequence 
listed in an above section. I'his section is a portion of the output if reipiestcd hy the user at 
input. 

.V.S.2 Sample Cases 

I wo separate cases are presented in the sample printout to illustrate the option ol Whether 
the microphone coordinates are cahicated hy ma.xtmi ing the determinant of the coherent 
\\.;\e equation system m whether th. \ are speufied ■ ither arhitrarily or as a set of existing 
microphone locations. Ihese sample cases demonstrate the execution of each option with 
d.it.i listed in l igiire (t. I hc length units in the printout are in centimeters; the time units m 
seconds. 
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1 lu' ^ll^l N.tmi'U' vMv' ilhiNti.ili'N dll’ o|'ii»'n I't ili'li immiiit’ imi u>|'lumi’ UuMlion^ loi .1 mIii.i 
Ill'll Ill'll’ lliii'i’ lUi'ili'N ail' I'li'i’a'Ml nil’ ill .1 lull iiii'li'i liiamoU'i aiiiiiilai iliul. Mu ti'i'lu'iu’s 
ai itii'M’ li'i'alii'UM'aii I'o um‘i 1 lo I'lMain ai'iuislir i'U'-imiiv iiu'.i'iiiu'iiu'iiIn. at a lu'iiiu’iu'V ol 
('.'Oil Ill’ll,', li'i lU’li’iiiiiiiiii!*. llu’ ampliliiili' aiul pli.iM' ol llu’ ( I.Oi. 1 l.l I. aiui i l. 'i iiiinli"i. 

Ml i ii'i'lu'iii’ I'l.ii'i’iiii'iil ll1ll^l I'l' ii’si t u'li’il ii' wilhiii llu’ j'.i'i'iiii'ii'ii’ I’l'iislraiiil s ol llu’ U'nI 
lai ililN . Ill tills Mliiation llio iiiuToi'lioiU' lOOtiliiiaU's aii’ I'oiiiuU’ii williiii a .Ul-i’iu axial aiul 
a .'10 lii’t'U'i’ I’lii’iimU’iviili.il M’aii’li latii'i'. Oiili'iil liom llu’ Mii ioplioiu’ 1 oialioii I’lopiaiii 
Tor this s.iiiipli' I’asi' I’ompriM’s llu’ sloi liaslii’ siuii li I'aK iilatu'iis. ilu’ itipiil I'aiaiiiiMiTs. aiul 
till’ I’.ili’iilali’il iiiii’ioi'lu'iu' loi’alioiis. 

I lk' slov’liaslii’ siMi’i'li pioi’i'iliiii’ lias l Otiu'Vjkkl lo .1 iiiK ioplioiii' li'i’alioii wlu’ii llu* I’liaiipi' 111 
Ilk' ili’li’i iniii.inl I'l'i’oiik’s as small as possibK’ lor siii.ill I’hani'.i’s in iiik ioplioik' loi alioiis 01 
w Ill’ll Ilk’ tnaxiimiin iiiiml'i’f ol soari’lk’s lias I'l’i’ii loiiii'li’loil. In this lasi' holli iiiii’ioplioiu’ 
loiMlioiis ill’ll’ riiiiiu'i! In tlii' pioii’iliiii’ lOini'n'i’il in Ino soaivlk’s Upon I’oiivi’nu’nii'. tlu’ 
llitily lai’s’.i’st v.ilik's of till’ ili’li’iniiiiaiit am! Ilii’ loiii’spoiiiliii!’, iiiii’iopltoiu’ looniinati’s laxial. 
I'ii’i’iimri’ii’iitial. ami lailial) aiv lisU’ii. .Mso. Ilk’ I’oiklitiotiiiij’ mitnhi’r is lislml loi liii’ lolii’r- 
i'iiln\a\i’ I’liiialioii sxsli’in ili’M iil'i'il l'\ Ilk' l Uiiiiilalivi’ si’t of niii’ioi’lioiU’ loialions. A I’oii- 
iliiiotiiiii’. luiiiilvr of .U 5 was ol'laiik’il for ilk’ lliri’i- wall iiioiinli'il iiinioplioik’s. imliialinii 
llial till' ili'li'niiinaiil of tlii’ I’liuatioii sysli’iii is woll away from /I’lo. lluis. Ilii’si’ lliivi’ loi’a^ 
li'.'iisaii’ ai'i'i'ptal'li' foi I’ali’iilaliii!; tlii’ nioilal slniiUiii’ by inalrix itni'isioii li’vliiiiiiik’s. 

\ lisliiiii I'f llii’ iiipiil variabli’s fi'llowi’ii lliis porlion ol tlu’ I'lilpiil. I liis lislitii’. iik’Uiiii's U’st 
ik’otiii’l I A ami I’oiiililk'ii paranu'li’rs. sloi’liaslii’ soaii’li I'ar.inii'li’is, aiul iiioiU' liunip. 1 ho 
iiitoinii’iiiali’ raliios assooialoil with llio iiioilo s'loup aro also labiilatoi!. I lioso valiios iiiolmli’ 
Ilk' axial wa\o iiiimbor in imils of iloiMvos-por li'iij’.llu llio oiooii \aliio. lx n,,, aiki llio \aliio 
of till’ oiiu’ii loi’tor 1 at tho iiiioioplioiK’ looalioiis. 

1 ho I’ooiiliiialos of tho solooloil mioroi'hono looalioiis aio lisloil. I hoso iiik iophoiios woro 
plaoi’il on llio iliii’l wall I’l'rrospoiiilitis', lo a railial oooriliiialos ol ’x 0111. 1 ho lirsi looation 
was spoi’ifioil at a ii'loroiioo looalioii willi /oro 0111 axial ami . 010 iloi’ioos oiii’iiinloronlial 00- 
oiilmalos. I wo looalioiis iilonlifioil b> tho I'toiitaiii aro plaooil al llio saiiio oiroiiniloioiilial 
liualioti as tho I'lisl tiiiorophono. I ho plaoi'iiiont of iiiioioi'lionos for iiioilos w ith Iho samo 
oil I’liiiiforoiil iai I'lilor is roslrioloi! Ii' iho satiio oiroiiiiili’ioiitial oi'orilinato booaiiso Iho iiiaxi- 
mi/ation of llio ilolormiiianl assooialoil with Iho oipiatioii s\slom is iiiiloiH’mlonl ol this 00- 
oiiliiialo. I ho .ixial looalioiis of Iho two solootoil niioiophoiio looalioiis aro ,U 1 0111 ami 1 > 0111. 

Iho si’ooml saiiiplo oaso ilhislralos tho option of spooifs iiij’. Iho iiiioiophono looalioiis. llio 
oaso was oxooiiloil with lliioo flush iiiountoii wall iniorophom’s al Iho saino oiroiiiiiloronlial 
oooiilmalo. I ho axial ilisl.moo botwoon oaoh imoti'phono was s|voifioil lo bo ton ooiilitiiolois. 
I his oaso is similar lo llio fiisl oaso boo.iiisi’ iho saino iiioilo oiinip, ltoi|iioiio\ . ami !’.ooniottio 
oonfiitiiratioii w oro usoil. 

Ihilpiit from Iho Miotophotio I ooalion I’loi’iam for this si-iomi saniplo oaso ooiiiprisos Iho 
ooiiiiil ioiiiii!’. iiiiiiibor ami oipoii \ahios of tho oi|iialioii s\ sioni assooialoil w ilh Iho siiooossno 
iiilioiliiolion of oaoh mioiophono looalioii. \ t mal oomliliomii;’ iiuinboi of ' (> w as obtainoil 
for tho lliroo iiik lot'hoiio looalioiis. iiulioalin,!', thal ihoso looatiotis aio .ilso aoooplablo loi 
ilolormiiiiiis* Iho aiiiphliiilo .iiui ph.isos ol iho i -l.l'l, i - 1 . 1 ). i I. 'I iiioilo sii lu lino 1 lio iiipul 
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H.-., "iIua'i a m.mlvr l.u 

aAi,.,.na.a, 

itl ^Miulul.iU’ K'l.ilions ..m ^ j | uci'pl;il'K‘ U'l l'‘H' ‘"'Utul mmU' 
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xiiikiinv ai'lv'iimiulinn. 


.V(, M.U lllNl. U! QimU MIM S 
llu' Murnpl"''''’ ^ 
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APPENDIX A 


Calculation of tlie Characteristic Numbers 


The cliaractcristic numbers K ® and (j£,, are defined to be the paired roots of the simultaneous 

m u 

equations 


< a 


[-4- j^CKm^r') + r')] »=i 

dr mju dr * 


( 


, a 
'm/i 


= 0 


O) 


r'^ + Q^" -p Y^(ok" r')lr'= 1 =0 

dr - V* m/Li mju dr mju * ‘ 


(2) 


b'or a given circumferential mode number, m, radial order, m. and hub/tip ratio, o, (where 
0 is not equal to zero); and Yj^ are the Bessel functions of the first and second kinds of 
order m. 

The following relations are used in the formulation of a solution 


d ' dx 

Jm (x) = V (x)^ 


(3) 


d ' dx 

*^Y^(x) = Y^(x)^ 


(4) 


2m 

■^m+1 ■■ "^m-l 


( 5 ) 


/ 1 

Jm ix) ■^m+1^’*^^ 


( 6 ) 


2 x "^ni '^m*l 


m 


= Vl 


puvkjeding page be 


ank mot htu^ 


25 


1 


(7) 


y'„(x) - I j;,, (X) .'' '11 ‘A! 

irx.l|„(\) J|i, 0^1 


/( s .l|„ (X ) 


Ml 


"1 . 
(X) - .1 


III (x)l 


, (I 0 

I I'tliii}’ K k :iii(l (.) - 1.) , Mini I'viihiiiliim ;il r 

lll/i IIIJLI 


Im 

J|l, (x) 

-■ I ; ( n ililtl ( 2 ) boaimo 


J' (k)k+(.)Y' (k)k-O 

III 

( 8 ) 

(o k)ok + Q Yj„ (ok) ok ^ 0 

t‘)) 

•I'mCklk 

1 rom ( 8 ), Q - suhstiUiting into ( 9 ) yieUls 

1 (1^ (K) K 



x ^ ' - !v! - r k -0 ( 10 ) 

^ m i ^ ^ 

I,ctr(K) = j;„(ok) y;,|K)o K“ j;^, (K)y|„ (nk)alf“ = 0 ( 11 ) 

Usini* tho expressions in (5). (0), (7). anil (11) then: 


Y 


f(K) = aK^IV,0,K) ;;^J„,(«K)| | r | ^ 




'* m )*' lx) / 

I 0(1 


111 


l(K)- a k 


>1 i - Um -1 (ok)- Jin ■'■'I ''Jni-I<*^> J,„(K)| 


n k J,„ (k) 


ir 0 K .1,,, (o k ) 




(M) 


III 


|:.,„.HMnll.n,s.v.,lu;.U^I^olv.,lm•solK, I. X ■ - . mUil nKj) KKj I)- 0 

loi sonu> 1 . A iMoaulmv .■nipKn mr, a o>ml.inal.oii ol' lalso pos.iion ait.l hiMMioii ir, hnu|.u-s 

is tluMi iisi'il li> ol'lain a viiliiv ol I 1* 

llavmc. L-alnilalo.l a v.iliu* ol K K/,;,,. llu' .•onvspo.ul.ni-, vaUu- olA) . an In- .alvnlalod. 


avni|’. I'.iH ii.iinn ■> ''inp 

Conil'inin)'. (S> aiul CM yii’lils. 


u' IK) ( .i;„loK>o 1 K . 0|V;„(K) » V',„O.Klol K 0 


iin 


from wliioli 

0 


J,„fK) I J„, foK)n 

v;„tKn V;,,(nKfo 


flM 


l-oi o 0 O 0 aiul K “ 0 is .li-fimul to K‘ ll.v rool of 

' ni/i 


It-, '^’ 1.' I ‘‘ 


Mr' 


n.fi 


I l(M 


li‘tlinp, K K '* . anil ovalualint’. at i' 1 . 1 1 (A Ivion.i's 


mp 


If UK) .I,,,' iK) K 0, tlu'i. fotyivlils 


im 


m 


UK) 


I IS) 


nation US) is oval..ati-ii for valiu's of Kj m t -Hi I ) ; . I . .1. • ■ • >"<!il f iK j) f (Kj j ) ^ . 0 

for sonu- valiu’ of i. A pioi ialmv mnploymi’. a amil'inalion ol lalso position and bisivlion 


toi-lumiiu's is tluMi iisi'd to obtain a vahu- of in llu' interval tKj |. Kj). 


APPENDIX B 


MICROPHONE LOCATION PROGRAM 
PROGRAM LISTING 




PAGE BLANK NOT FIT.MEO 



OS' ponir . 
^^Oli qVAlAT 


PRAT1 C HHllNer AIRCRAFT DIVISION 

sc.fanlid.la 


VER 

9.0 


07/2S/77 

i 2 .:>o.oo 


♦WRITE PRINT, T99«#01 

DATA SIT T69901 AT LEVEL Oil AS OF 

DATA SLT TA9901 A1 LEVEL 008 AS OF 


06/16/77 

06/06/77 


00001 
00002 

THIS PRUU^AH CAICULATES THE ‘BEST* PLACEMENT FDR HICROPHONES IN A OUCT00003 

00006 

EIOVECISO, 501*00006 
00006 


DIHENSIIIN AA(200), INOIIO), WQRKI5IO0I* E1GVALI50)* 

I 10(50,50), AMPL(50), PHASEI50) 

I XTEKNAL HI 1CAL 

Cl-MnON /HAIRX/ NUIHI, ONI, RO, DET, HATRIXI 50,50) 

CCM.U.N /I'kUblS/ IWALL, ISAHE, NPKUbE, KNPROB 
ComMCH /CNSTNT/ NMIKiS, NHOUtS, SILMA, 8, MX, FRLQ, A, OHEliA 
COM;>uH /IMUS/ fcHU(50,50), ItHpRT 
CLHMi.N /Mimi/ KHO)50), 0MU(5O) 

fj /hl.l IS/ MODI 150), MU(50), IWAVE(50) 

COPlUiN /HA VINO/ K.XI50) 

COMM./j /rtIKtS/ Nl)JM, XI50), k(50), THETA(50) 

U:.n oN /Si ARCH/ hVAR, IPKNT, NSRCH, NLVALt lEXPNT, NTABLE, NSTAT 
tlif.UifJ /bLlSL/ UUMK5), PI 

common /l.UUNuS/ Xb(jUNU|2,50), RBUUND (2 ,50 ) , TBDUND (2, 50) * XA, XB* 
' TimlN, TKMAX, XLIH, RLIH, THLIH 

(>I100CN 

OI-I.KAlt, KaDDIG 


/UINOTN/ 
/ANCI IS/ 


C 

c 

c 


c 

c 

c 

c 

c 

c 


COMMON 
COHMuN 
RIAL KMU, MX 
COMRUX UHL, ZERO, 
1 tlGVEC, LO 


DET, MATRIX, FACTR, OIVSR, EXPNT, KX, EIGVAL, 


INPUT CASt VAKlAbLbS 

20 CALL INPUT) UNO ) 
IFI lENU .GT. 0 ) 


GO TO 9999 


CALCULATE CHARACTERISTIC NUMBERS KHU AND QMU FOR EACH SET OF 
CIRLUMI tHLNTlAL MODE NUMBER ANU RADIAL ORDER 

CALL KGCAL 

CALCULATE AXIAL HAVE NUMBER 


FILH B OMEGA / A 
AMACH * I. - MX * 
liU <>U I-1,NMUULS 


MX 


RAlJlCL 
IH KADICL 
25 KXd ) 

1 

30 Axn ) 

1 

AU CUNTINUE 


) 


- FLOW 2 - 

ChPLXI -MX 
‘.WHTI ABSI 


AMACH 


/ B I 2 


♦ FLOW 
KADICL 


CMPLXt 

AMACH, 


( -MX 

0.0 ) 


* ( KMUIl) 

25, 30, 30 

/ AMACH, IMAVE(l) • 

) ) / AMACH I 
GU TO 60 

MUH ♦ IHAVMI) • SORTI RADICL 1 ) 


OITLKMINL MICRDPHUNE LOCATION COMPONENTS 
SIOUiASnc SIAKCH TLCHNlUUL 


, X, R, AND THETA, USING A 


00007 

OOOOB 

00009 

00010 
OOOTl 
00012 
0G0I3 
00016 
00015 
000 1 b 
00017 
OOOIB 

00019 

00020 
00021 
00022 
000 23 
00026 

00025 

00026 
00027 
0002B 
OOOrt 

00030 

00031 

00032 

00033 
00036 

00035 

00036 

00037 
00U3U 
00039 
0006U 
0001,1 
00062 

-00063 

00066 

00065 

/00066 

0UU67 

OOU6H 

00069 

00050 

00051 


preceding page BTiANK NOT FITjMED 


n 


I 


ooo or", oon 


PPATI L HHIINEV AlRCRAfT DIVISION 
&C.PANL1H.L4 


VER 07/?5/77 

9.0 12.t>0.00 


120 IHALL 
IPKUbk 
HPRUBt 


> 0 
= 0 
» 0 


Sti UP FlRSI ROW OF HAIRIX ANO INITIAL DETEftHINAMT VALUE 

CALCUlATt CHAHACURISTIC E-FUNCTION VALUES ASSOCIATED WITH R 

HPRIHt « Rill / B 

CALL tHUCAH kPRlMlt FMUIlfU* NHODES > 

UU lAO J«ltNHUUES 


LXPNl 

1 

HATKlXU.Jl 

i 

lUI 1 iJl 
lAO CONTINUE 
btl 


CMPLX< O.Ot REAL! RXUl I * Kill ♦ HODEUI * 
ThCTA(l) ) . . 

^ EMUJJ.U • CLXPI EXPNT ) • £XP| -AIMA&I KXIJl I 

xm I 

= MAlKlXl ItJt 


c HAiKixn.n 

SET UP ARRAYS FOR STOCHASTIC SEARCH 


INDt II 

INU<2) 

1 N 0 I 3 I 

1N0< 

iNOtSI 

1N0(6| 

INiitVl 
INM 10 1 
AVO 

MKSTKT 


NVAK 

NSKCH 

NkVAL 

N1ABLE 

NSTAT 

0 

0 

ttXPNT 

IPRNT 

I 

.9 » I XB - XA } 
NUIH * 1 


DO 440 JJcHKSTRT.NHlKES 
HKII E(6.9000l 

9000 FOKHAl (IHl, T44, •*** MICROPHONE LOCATION COMPUTER PR06RAM 
JJ .Lk. NMlKkS I CO TO 1*7 


C ALL MlCROPHUNk LOCATIONS HAVE BEEN INPUT. CALCULATE CONDITION NUMBER 
C 

DU 430 JJJc1,NH1KES 
J - JJJ 

CO TO 17 5 

IK iPRNT »GI. 0 I WRITE (6.90011 4 

9001 FORMAT ( /, IX, •... STOCHASTIC SEARCH OUTPUT FOR MICROPHONE », I2I00101 

f 0 0 1 02 

C IF MODI NUMBERS CHANGE, RESET ISAME INDICATOR ANO ATTEMPT TO PLACE 
C HIKE uN DUCT HALL 


00052 

00053 

00054 

00055 

00056 

00057 
00050 
000 59 
00060 
00061 
00062 

00063 

00064 

00065 
• 00066 

00067 

0006B 

00069 

00070 

00071 
0C072 

00073 

00074 

00075 

00076 

00077 
0007U 
00079 
QUObO 
OOOHl 
00002 
000B3 
00064*11 
00085 

00066 
Q00B7 
UOOHB 
00089 

**••00090 

00091 

00092 

00093 

00094 
0009 5 

00096 

00097 
0009B 
00099 
OOlUO 


00103 

00104 




n r> r\ 


page £9 

Poor quautyj 


PKA11 C HHlTNey AIRCRAFT DIVISION 
SC.RANLlb.L'* 


VER 07/25/77 

9*0 12.50.00 


1H( ISAHE .tU. 0 ) 

1 E( HOUMJI .tU. HUOEtl> ) 


UAHE 
1MAI.L 
150 r^Lilhl 
N(ilH 

155 AA(1 ) 


« 0 
» 0 

c J - I 

= J 
» XA 


1E< IHAIL .Nr.. O 1 
AA <2 » * TMHtN- 


IM SILRA t 
160 AAlil 

165 AA(3 I 
170 AAlA ) 

AA(5 ) 

AA16 i 


- .1 

« B * SIGHA 

c Xb 
•= IHMAX 
•= B 


GO TO 150 
GO TO 150 


AAIU > XA * AVG 
160. 165 


160 
GO 10 170 


SET UP bliUND ARRAYS 

Xbl.tlOU(l»J) - AAI l» 
xiii.uriu • j I aA(a) 

TbCjUNUlI.J) AA(^) 

TlCJIJOD(/,J| .. AA(5l 
R(:(ajM)<l,J) = AA4 3I 

« AA16J 

CAti. sriK.Mi DPtLALi INO. AA I 
SIT UP HiCKOPhUNt LDCAllUN COMPONENT ARRAYS 


X(J> 

IhkI A( Jl 
k(J) 


AAUt 
- AAJ31 
» AA<3| 


SET UP HUH J UP MATRIX USING THE VALUES UP X* THETA. AND R 

CALCULATE CHAHACTERlSTlL E-FUNCTION VALUES ASSOCIATED WITH R 

175 KPKIHE • KIJI / U 

CALL EHUCALI HPRIMLi EHUU.Jl* NHOUES I 
DU 150 1*^]|NM(1UES 

tXPNT B CMPLXI O.d. REaLI KXI 1 1 I • XU) *■ MODEUI * 

1 TUFTAIJI 1 

HA1KUU.il = Emin.Jl • CtXPI EXPNT I * EXPI -AlMAGl KXIll 
1 XUI I 


EUlJit I a HAThlXIJ.n 

150 CCNTINUL 

1F( Jj ,GT. NHIkiS I 

1R1ANUIL AK12E THE J a J HAIRIX 

FALTR - UNE 

DU 220 Ul.NUlHl 


CO 10 230 


00105 

0010b 

00107 

OOlCB 

UOIOV 

00110 

0011 1 

G0U2 

00113*11 

DOHA 

00115 

00116 
00117 
OOllB 
00119 
G012U 
00121 
00122 
00123 
0012A 
OU125 
00126 
00127 
0012'i 

00129 

00130 

00131 
GUI 32 

00133 

00134 
00136 

00136 

00137 
0013B 
00139 
OOIAO 
00191 
0U192 
00193 
00199 

00195 

00196 

00197 
* 0019b 

00199 

00150 

00151 

00152 

00153 
U0159 
0015 5 

00156 

00157 




* ‘ • CRNJGI HATK1XIJ,|) ) / 

IlJVSR - „ *1*’ •••?»• fAUTR ' 

^ CA8M JJiJjJ HATRIXIJ,!, , , 

iSTKT ^l4-i 

00 200 Ki'l&TRTtNRdUtS 

200 LOmiNUt**^* * «*TRIXI1,K» ♦ DIVSR • HATKIXIJ.RI 

^ 220 COfiTlNUt 

C 

c CALCUUTl CONOItJON NUMBER Of THE MATRIX 
230 MIRENO = j 

^ ■•‘IL ur.cc, tc. ^ 

eigm]* “ eicvAuii I 

n,, “ EIGMIN 

T?7^^ ® CABS! EIGVALdI I 

EH EIGEN - EIGHIN t 
2A0 ElGHlN c eigen 


OIVSR 


260 IM Eir,tN _ eigmax » 

290 EIGMAX t EIGIN 

300 CUNHNUt 


2A0t 260 t 260 

00 TO 300 
300, 300, ?B0 


CONONO 


= EIGMAX / EIGMIN 


07/25/77 PAGE SERI At 


12.50.00 

O015H 

00159 

00160 
00161 
0016? 
00163 
00169 

00165 

00166 
00167 
0016B 

00169 

00170 

00171 

00172 

00173 
00179 

00175 

00176 

00177 
00179 
00179 
OOlBO 
0018 1 
00182 
00183 
00199 

00185 

00186 
00187 


. 'HASA 0. A.«K,AA..., AA.X, 

00189 

PMASeO J “ EIGVALCLI 1 00190 

310‘goNhI ' B 

*. .2, . oSiH 


2'amhliiuue VnisE. 

DU 315 I.H,M1KENU 

W03 rOR/!AH’Al!r*ip’ ‘^HASUU OOIV^ 

c 1 r,„ ,0 ooloy 

c IE CUNU.riUN number ,S .onou., CUNIINUE. OTHERH|SE PUT MIKE ON PROBE SS 

^ 320 !f! SrrT!T5" ’ --ulSi -o? 

c hike is EITHER UN HALE UR NU| THE MHST HIKE UN A PROBE. IN AN7 EVENT.^Hf 


OiRlGINAi; PAGEl la 

Otf' P(K)R QUALmH 


r 


viKH'i'isAi: r/uii’. t;. 
' >1- ujoit 


PRATT fc WHllNkV AIRCRAFT DIVISION 
SC.PaNLIB.L^ 


VER 


07/26/77 

12.60.00 


C TRY AND PUT MIKE ON ANOTHER PROBE 
C 

IPRUBfc » l ♦ I PROBE 

^ 1F< IPROBE .LE. NPROBL ) GO TO 3A0 

C tRKUll - iHt NUMBtR OF PROBES EXCEEDS THE MAXIMUM ALLOWABLE 

3-TO WRnEI6,1000) 

lono MIKHAII //, IX» the number of probes REOUIREO EXCEEDS THE 0021V 

IMAAlMuM KlRHllTtl). EXECUTION WILL- BE TthMlNAlED } 00220 

, GO TO 9VV9 0022 1 

SaO MPHUl t «. 1 

IHAl L B -1 00« 

IWAl.L L. 1 0022 3 

^ GO TO 155 0022A 

r ACCEPTABLE. IF ON HALL, CONTINUE. IF UN PROBE, READJUST PROBE 00226 


00211 
00212 
00213 
0021 A 
00216 
00216 

00217 

00218 


C 

c 


1 NFLKMATIUN 

360 1F< IHALL .EU. 0 1 

HFKUBi « 1 ♦ MPROBE 

1F< HIRUBE ,GT, MNPROB I 
iHALL e I 

380 IPRIlbl * I «• IPROBE 

1E( IPRUBE .GT. NPROBE ) 
IWALL K -E 

HPkUBE » ] 


GO TO A20 
GO 10 380 
GO TO A20 
GO TO 330 


C CALCULATE The UETERMINANT of THE J X J MATR 
C 

OET B MATRIXIJ,JI * PACTR ♦ DCT 

GO TO VAO 

A30 CONTINUE i 

AAO LTNTINUE 
C 

C PRINT RESULTS 
C 

CALL »‘R1NT 
C 

C RECYCLE FOR NEXT CASE 
C 

^ GO TO 20 

9VVV STOP 

END 

^ SUhkUUTlNE ANGPOST ANGLE, NUMBER I 

C THIS SUUHOUTINL CONVERTS NEGATIVE ANCLES TO CORRESPONDING POSITIVE 
C ANGi. I S 
C 

OTMENSIUN ANGLEO) 

DATA UlGREE / 360.0 / 

C 

1)0 80 UliNUHBER 


0022 7 
00228 
Q022V 
U02TO 

00231 

00232 

00233 
0023A 

00235 

00236 
002 3 7 

00238 

00239 
002A0 
002A1 
002A2 
002A3 
002AA 
002A5 
002A6 
002A7 
002A8 
0o2A9 
002 50 

00251 

00252 

00253 
0V25A 

00255 

00256 
00^51 

00258 

00259 
0026U 
00261 
00262 
00263 


i-. 


1. 


o r n 


PRATI & MHllNLY AIRCRAH OIVISION 
S1..PANL1H.L^ 


VtR Ql/Z^m 

9.0 12.^0.00 



IM ANGLE 

( 1 1 

I 


20 

OU AO J~1 

,10 




OILT A 


s 

J » DEGHLE 


IM ANGLL 

1 1 1 

♦ 

DELTA 1 

40 

CUNI INUt 




60 

ANt.l MU 


s 

DELTA ♦ ANGLEll I 

no 

CUM INIIL 




999 9 

Kt TlIKN 





1 NO 





FL.NC nUN 

BtSL 

11 

X 1 


20, 

eo« 

ao 

0026A 




00265 




00266 

AO, 

6 0 1 

60 

0026 7 


002 6 >i 
002A9 
00?70 
00271 
002 72 
OOi’ 1 3 
002 7<i 
002/5 


THIS KINlIIUN CAtUILATES VALUES OF THE EOUATION DEFINING THE 
DlFURtNUAL lOUATIONS FUR A N0N-2tR0 HUB/TIP RATIO 


LUHHIW /IIISSL/ ISION, JSIGN, OELKMU, TOL, 
LuHMl.N /LNSINt/ 0UM1I21, SIGHA, 0UH2151 


Ht PI 


XI 

CAll 

CALL 

CALL 

CALL 

CALL 

CALL 


OLSJl 
BLSJl 
HISJ( 
bt SJI 
Bl SY( 
1 .lSY< 


= X • SIGMA 

Xlt M-.ISIGN, IMJHLt TUL. lERl 
XI» H . LMjXli TULt UH2 ) 

X« K| LMJf IlILf 11H3 I 
X, M-JSIGN, LMJPl, TOL, ILRA 1 
Xt Ml LMYXi 1LR5 » 

XL, H» IMYXlt 1LR6 1 


LMJMl 

iHjXl 

LHJ 

IHJPI 

IHYX 

LMYXI 

A1 

A2 

A3 

AA 

A5 


JSIGN 
1 S I I.N 
IStl-N 
JStl.N 
1 SI I.N 
ISIt.N 


LMJMl - 
LMJPl - 


2. 

2 . 

Al 


A1 

A2 

A2 


ISIGN * LMJMl 

I HJxl 

LMJ 

ISIi.N * EHJPI 
LMYX 
I MVXl 

( M * JSIGN / XI 1 • EHJXl 
{ H * JSIGN /XI* EHJ 
/ ( PI * X * EMJ I 
/ ( PI ♦ XI ♦ EMJXl I 
• ( LMYX / LHJ - LMYXI / EHJXl I 


BLSLl 
KLIUhN 
LNl) 

LUNCH UN bLSL2l X 


XI • X • 1 A3 - AA ♦ A5 I 


I 


PI 


C THIS HJNLllON CALCULAILS VALUES UF THE EOUATIUN DEFINING THE 
C imHKlNllAL LWllAlIUNS FUR A MU0/11P RATIO OF URO 
C 

LUMMUN /ttlSSL/ ISU.N, JSIGN, OtlKMU, TOL, 

U.MHON /LNSlNI/ UUHII2), SIGMA, IHIM2rj) 

CAIL BISJl X, M-JSIGN, LMJMl, UlL , ILRl I 
CAM lUSJl X, H, LHJ, TUL, lER2 I 


SYSTEM 01 002 Hi 
00277 
002 /H 
002 79 
002110 
U02H 1 
002 h 2 
002h3 
002BA 
002bS 
002H6 
002B7 
002 (IB 
002(19 
002VO 
00.91 

00292 

00293 
0029A 
0029S 

* 00296 

00297 

0029 B 
00299 
OOJOO 

0030 1 
003u2 

0030 3 
OOUlA 
00 30 b 

00306 

00307 
system OFGOjCH 

U0jU9 
UU 3 1 U 

0031 1 
00312 
0031 3 
003I-, 

00315 

00316 


c p r> n r r> o r>nrirn r> o r 


I 
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PRAIT I WHlTNty A|RCR*fT DIVISION yER 07/25/77 

SC.PANLIB.LA »,0 12.50.00 

CNJHl • JSIgN tSlCN * EMJHl 00317 

EMJ B I&ION « ENJ 0031H 

BESL2 X • EMJMl > H • JSICN * EHJ 00319 

RtTlIKN 00320 

INU 00321 

SUbklHillNt BEiJl Xt Nt BJ, 0» lER I 00322 

t 00323 

C THIS SUItKOUTlNE r.ALCULAlES ThE J BESSEL EUNCTION FOR A GIVEN ARGUMEN1t0032A 
C X> ANii UKUkH N. THIS SUbKUUTlNE MAS-TAKEN FROM THE IBM SCIENTIFIC 00325 

C SUhRUl/llNfc PACKAGE 00326 

^ 00327 

bJ 0.0 00328 

Ikl N *Gk. 01 GO TO 20 00329 

00330 

C ERROR - NEGATIVE ORDER. SET ERROR INDICATOR TO 1 AND RETURN 00331 

t 00332 

ItR = 1 00333 

GO TO 9999 00339 

20 1F( X t 30» 60 00335 

30 1F< N .GV. 0 ) CO 10 90 00336 

bJ “ 1.0 00331 

GO TO 9999 00336 

00339 

ERROR * ARG'JMfcNT ZERO OF NEGATIVE. SET ERROR INDICATOR TO 2 AND R^TURNOOSAO 

00391 

AO lER » 2 00392 

GO TO 9999 00393 

00399 

CAICULAIE MAXtHUH ORDER NUMBER THAT CAN BE PROCESSED FOR X. 00395 

IE X .LE. 15. N RUST BE LESS THAN 20 * 109X - X9«(2/3I 00396 

IF X .GT. 15. N must be LESS THAN 90 * X/2 00397 

00399 

60 IFI X - 15. ) 00. 80. 100 * 00399 

60 N1ES1 * 20. ♦ 10. » X X •* 2 / 3. 00350 

GO TO 120 00351 

100 NIEST e 90. « X / 2. 00352 

120 IF! N .LT. NTEST ) GO TO 190 00353 

00359 

ERHOR “ ORDER RANGE COMPARED TO X IS NOT CORRECT. SET ERROR INDICATOR 00355 
1U 9 AND RETURN. 00356 

00357 

lER « 9 00358 

GO TO 9999 00359 

l90 ILK & 0 00360 

Nl » N ♦ I 00361 

8PREV ■ 0.0 00362 

00363 

UmPUTE STARTING VALUE OF M 00369 

00365 

1F( X - 5. I IbOt IBOt 190 00366 

160 MA K X ♦ 6. 00367 

GO TO 200 0036S 

180 HA ■ 1.9 * X F 60. / X 00369 


ry Cf ^ frtoo n n ^ nr*c^ 
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PKaU L HHllNJt AlKCKAfl OlVISlUN VtR 07/P6/T? 

SC.t' ANU 9*0 12«!»0«00 

?00 MB R N ♦ IHX( X ) / 4 * 2 003YU 

H7LKU *- MAX0( MAi HB ) 00371 

00372 

iL I llKJ-l R I inn Oe M 00373 

00374 

MHAX « NJbST 00371> 

220 01) 320 Hi.n/tK0,MMAX,3 00376 

003 77 

SI T HM» , HM-l I 003/a 

00379 

bMl t. l.0k*-26 003B0 

IM 0,0 0030i 

ALHIIA » O.O 00332 

i 00jh3 

IM ( H / 2 > * 2 ,F.U* H ) J1 k >1 00364 

M2 « H ~ 2 003B6 

1.CJ 2B0 K-1.M2 0G386 

M^ = M - K 003ti7 

HMK - 2. ♦ FLOAT! HK » * FMl / X - FH 003bB 

FM a. FHl C03b9 

f Ml * OHK 00390 

IM MK - N - X > 260» 240» 260 00391 

240 l-.l BMK 00392 

260 JT = -JT 00:193 

S a 1 ♦ JT 00394 

alpha b alpha * 6HK * S 00396 

2a0 (.t.NTlNUt 00 396 

00397 

HMA b 2. • FHl / X ~ FH 00398 

IH N .LU. 01 bJ « bHK 00399 

alpha » ALPHA ♦ BMK 00400 

8J = bJ / ALPHA 00401 

in ABSI HJ - HPhtV I - ABM D * BJ » » 9999« 9999t 300 ' 00402 

300 bPKbV B BJ 00403 

320 Lt.NTlNUt 00404 

00406 

LRKIJK - ktUllIKtO lULLHANCE NOT UBTAlNtO. Stl ERROR INOILATOR TO 3 AND 00406 
kl I liPN 00407 

0040B 

ILK “3 00409 

9*799 KtlUKN 00410 

LH() 00411 

bblthUUriNt BF$Y( Xi Nt RYt ILR I 00412 

00413 

THIS SUBROUTINE CALCULATES THE Y BESSEL FUNCTION FOR A CIVEM ARGUMENT ,004 14 
X, AflU LKUl R N. THIS subroutine haS TAKEN FROP THb IBH SCIENTIFIC 00416 
SOBRuniNI package , 00416 

. 004 IT 

Ilk •• 0 0041 B 

in N ,GL. 0 1 CO TO 20 00419 

00420 

EkkliR “ NtGATlVE UROER. SET ERROR INUICATOR TO 1 AND RETURN 00421 

00422 


LS 


a ^ 
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♦•RATT t WHlTNfeV AIRCRAFT DIVISION 
SC.PANCIR.LA 


VER 

9.0 


07/?S/77 

l?.SO.OO 


ICK 

iO 1F< X > 


AO IkK 


CO TO 9999 
A0» AOi 60 


CO TO 9999 


BRANCH IF X IS I ESS THAN OR IQUAC TO A. 


60 IM X - A. > 
CALCULATt VO AND VI FOR 


80 


OOA23 
OOA2A 
00A26 
00A26 

ERROR ' ARCUMEMY ZERO OR NEGATIVE. SET ERROR INDICATOR TO 2 AND 

00Ai<9 
OCAAO 
00«3l 
004 3 2 
OOA33 
O0A3A 

004 j 6 
004 36 
00A37 
00A3B 
00A3V 

C I ( I -.00000370A3 • T2 ♦ .000017356S I ♦ T2 - OOAAO 
.0000487613 I • T2 ♦ .000173A3 I • T2 - .001763062 OOAAl 
» * T2 ♦ .3969A23 00AA2 

(III .U000032312 • T2 - .00001A207B » • T2 ♦ 
.00003A2A68 » • Y2 - .0000869791 I • T2 ♦ 

,u0u466A32A I • T? - .0E2A6694 

I I I I .UU0UU42AIA • T2 - ,OUOC20092 I • 12 ♦, 
.o00y^807J>9 ) • 12 - .000223203 » * T2 ♦ 

.00^^21826 ) • 12 ,3969423 

I I I I -.OO0UO36A94 ♦ 12 ♦ .00001622 1 • T2 - 
.0000398708 I • 12 ♦ .000106A7AI 1 • 12 - 
,00063904 » 4 12 . .03740084 
2. / S«h1 I X 1 
A * 11 

X - .7863982 , 

A • FU 4 SINI C 1 4 B • UU • cost C I 


T1 

12 

Po 


4 

Tl 


lOOt 100, 80 

X GREATER THAN A. 

X 
Tl 


1 

2 

( 

1 

2 

1 

2 


CO 


PI 


Cl 


A 

8 

C 

YO 

Y1 


-A 4 PI 4 CtlSI C 


1 * b 4 ui 
Cl) Tn 160 


4 SlNl C I 


CALCULATE VO AND VI FUR X LESS THAN OR EOUAL 10 A 


100 XX * 

X2 * 

1 

SOH “ 

TIRH “ 

YO * 

UU 120 L«l*16 
)F1 L .NE. I I 
I L * 

IS » 

UKH * 

1 

YO 4 

120 CONIINOE 

URN * 


.6 4 

XX 4 
ALOGI 
0.0 
I 


X 

XX 

XX 


I 4 ,6772167 


SUN 


1. / FLOAT 1 L~ll 4 SUM 


L 

I - SUM 
I -X2 • TERM 
IS I 1 
URM 4 VO 


/ I 


FL •• 2 I I ♦ I 1. ~ * 


- .% 1 


00AA3 

UOAaA 

004A6 

004A6 

00A47 

00440 

00449 

00460 

00461 

00462 

00463 
U0h6A 
00466 

00466 

00467 
004 6 6 
00469 
0046U 

00461 

00462 

00463 

00464 
00466 

00466 

00467 

09468 

00469 

00470 
0047 I 
U0472 
0047 3 

004 74 
00h7 6 


i>) 


XX • I T 
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SUH « O.U 

VI « URN 

nu 1^0 L>2«16 


VI 

lAO CONTINUE 
PI2 
Vti 

VI 


» 1. / FLOAT! L-I I ♦ SUN 

* I 

• FL - 1, 

‘ 1 - SUH 

“ I -X2 * ItRH / ( FL * rtl I ) ♦ I I TS - .6 / 
/ I Is ♦ .S / Fll ) ) 

= TtRM ♦ Yl 


00476 
004 V 7 
00471 
004/V 
004HU 
004H1 
004U2 
PL > OO^bi 


B *6366190 
“ P12 ♦ YO 

» Pli ♦ t Yl - 1, / X I 


C LHtCK IF ONLY YO OR Vl is DtSIRtD 

I6u IF| N *&T. 1 I 
C 

C RMURN YO OH Yl AS REOUlRtO 
t 

BY YO 

IF! N .EO, I I 


GO 10 180 


BY B Yl 
CO TO 9999 


C PERFORM KttURRINCfc llPfcKATIONS TO FIND VNIXJ 

180 VA s YO 

VB « Yl 

R * 1 

200 1 = FLOAT! 24R > / x 

VL = T * YB - YA 

IF! ABS! VC I - 1.0L70 I 240, 240, 220 » 

tXt-fcEDED HAGNlTUOt OF 10**70. SET ERROR INDICATOR 

220 ILR « i 

240 K . , , K 

IF! K .E«. N I GO 10 260 

VA B YB 

Yb ■ YC 

260 BY . YC 

9999 RETURN 
EEJd 

blUCK DATA 

UimhuN /OFAULT/ NIOC, HiR, OR, LMX, FRO, HISOK MUSI50), IDIR 
» r.ONDNU, ILNU, XMIN, XMAX, DbLxI DfcLTH. 

^ OELR, NO, XMISO), RM|SO», THMiSOI, |PT, NSRH, 

^ TMAX* NPkUO, NPROB, N1AB, NST, SPEtO, 


004b4 

004H6 

004U6 

00487 

0U4B8 

004B9 

00490 

004 91 

00492 

U049A 

00494 

00496 

00496 

00-.97 

0049B 

00499 

00600 

00601 

00602 

00603 

0U6O4 

00605 

00606 
00607 
006uB 
00609 

TO 3 00610 
00611 

00512 

00513 

00514 

00515 
00616 

00617 

00618 
00619 
00520 
00621 
0062 2 
00623 

ISO), 00624 
00626 
0062b 
TMIN, 00627 
006.; 1 


■r 


r 


1 
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PRATT C WHlTNtT AIRCRAFT DIVISION 
SC.PANL16.U 


i GtOMlTRY ANU-JtST CONDITION OEFAUCT VALUES 

: NLOC - NUHBLR OF HlCKaPHONe LOCATIONS OESIREO 

HIk - HUH / TIP RATIO 
: OK - OUltk RADIUS UF OUCT 

tMX - AXIAL HACH NUMUbR 

; FkU «■ TLiT FRtUUtNCY 

DATA NLUL / 2 /* HTK / O.A3B /, OR / S.O /, ENX / -O.OT /. 

I FRO / 3100.0 /» SPEED / 13566,24 / 

MODE DEFAULT VAIUIS 

M - CIRCUMFERENTIAL MODE NUMBER 

MUS - RADIAL ORDER 

lOlR - NAVC DlRLCllON INDICATOR 

DATA M / ?*-2, ab* 0 /, HUS / 0* 1« 4B60 /, lOlR / 2*1, Aa«0 / 

CONDITION NUMULR DEFAULT VALUE 

CONDNO - MAXIMUM CONDITION NUMBER 

DATA CONDNO / lUO.O / 

STOCHASTIC SEARCH PAKAMLTEK BOUND-DEFAULT VALUES 

XMIN - MINIMUM X VALUE OF OUCT 

XMAX - MAXIMUM X VALUE OF DUCT 

~ ****'■ ‘^*^TANCE ALLOWED BETWEEN MICROPHONES 

ri radial DISTANCE ALLOWED BETWEEN MICROPHONES 

THIN ~ H annular DISTANCE ALLOWED BETWEEN MICROPHONES' I 

THIN - minimum UOuND FOR TllETA 

TMAX - MAXIMUM BOUND FOR IMETA 

DATA / 0.0 /, XMAX / 12.0 /, DEU / l.O /, OELR / l.O /, 

1 OELTh / 11.5 /, ihIh / 0,0 /, TMAX / 360. / 

MlCROPHONt location DEFAULT VALUES 

ND - NUMBER OF INITIAL FIXED MICROPHONE LOCATIONS 
AM - AXIAL CDMPdhiNTS Of MlCK(JPhONF LOCATIONS 
RM - RADIAL CDHPONINIS IlF MlCHOPhUNE LOCATIONS 
ThM - AHLULAR CUMPUNtNTS DF MlCROPtlONL LOCATIONS 

^DATA ' S0*0.0 /, RM / 2*5.0, 48*0.0 /, 


stochastic seaklii default values 

IPT - print indicator fur SIARCHES 

NSPH ~ MAXIMUM NUMBER OF SEARCHES 


07/25/77 

12.50.00 

00529 

00530 

00531 

00532 

00533 

00534 

00535 

00536 

00537 
0053B 

00539 

00540 

00541 
0054 2 

00543 

00544 

00545 

00546 
’ 00547 

00548 

00549 

00550 

00551 
. 00562 

00553 

00554 

00555 

00556 

00557 

00558 

00559 

00560 
IOEG100561 

00562 

00563 

00564 

00565 
CC566 

00561 
00569 

00569 

00570 

00571 
C0572 
005V3 

00574 

00575 

00576 

00577 
CC578 
005/9 
005HU 
005B1 
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iC.I'ANLIH.I-A ».0 12.>0.00 

C NVAL - NUHBtH UF EVALUAT1UN& PFR SEARCH 005S2 

C IXPNT - EXPUNLN1 QF TULtRANtfc I E.G. TOL • 10. •• IXPNT » OOSB^ 

C HIAA > NtiHULK UF tUNLNIS IN UFTfcRMlNANl TAOLE OO^BA 

C NST - NUMhLH (IF UlNINTS 1u OF USED IN CALGULAIlNG S7A11STICS 00;>Ui> 

C NVAK > NUMUtK UF INDb I'L NUF N1 VAKlAbLkS 00;>b6 

C 005>fa7 

DATA IPT / 0 /♦ NSRIi / *>0 /, NVAl / *.00 /* IXPNT / -2 /, 

1 NTAB / JO /f Nil / JO / 00!>B9 

UlHHUN /ilAKLM/ NVAR, llUH(6| 00!>90 

Data nvah / j / oo«, 9 i 

C 00592 

C PROUL UIFAUALI VALUti 00593 

C 0059A 

C NPRUB - KAXlmjN NUMBER UF PROBES ALLUMEO IN OUCT 00595 

C HPKiiB - MaXIHUM number OF MICROPHONES PER PROBE ALLOMEU 00596 

C 00597 

DATA NPHUB / I /. MPR08 / 1 / 00598 

C 00599 

C GENERAL PKTNI Otf AULI VALUES 00600 

C 0060 L 

C I EMU - pRlNt INDICAIOR FOR CHARACTERISTIC E-FUNCTIQN VALUES 00G02 

C 00603 

OAIA lEMU / 0 / 0060A 

, 00605 

BESSEL FUNLIlUN VALUES 00606 

00607 

COMMON /BLSSL/ DUH212)t DELKHUt TQL« MM. PI 0060B 

OAIA liLLKHU / 3.0 /* IQL / .0001 /t Pi / 3.141593 / 00609 

00610 

ANGULAR LUNVEKilUN VALUES 00611 

00612 

UEGkAU ~ ULGKEES TU RAOIANS 00613 

RADUEG - RADIANS 10 ULliHEES * 00614 

00615 

CTHMUN /ANGLES/ UbGHAUf KAOOEC 00616 

DATA OlGKAl) / 0.D1/45JJ /, RADOEC / 57.29578 / 00617 

00618 

CliHPLtX ONE. ZERO. DUM3 00619 

CUHHUN /NaTRX/ NOImI, one. zero. 0UM3I 25011 00620 

DATA ONE / II. .0.1 /. ZERO / (0.|0.) / 00621 

END 00622 

SUitKUUIlNE IHUCALI RPRIHE, EMU. NMODES I 00623 

00624 

THIS SUBKUUllNL LALCULATlS NMUOES CHARACTERISTIC E^EUNCTIEM VALUES F0R0Q625 
A PAktl(.UlAK radial VALIIL, KpRlFl. 00626 

0062 7 

(ilHENSUiN LHUin 00628 

CGMMI.N /KONU/ KHUI50I. UMU|50t 00629 

UIHMU4 /NUDES/ NODE 150). OUN2|tUO) 00630 

( UHHUN /UtSSl/ ISlGNt JSlGNi DLLKHU. TOL* N. PI 006^1 

RIAL KMII 00632 

i.U 40 IkI.NHMUIS 006j3 

H ^ lAhSI H(rt>E(l> t 00634 


I 
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“RATI C HKIlNtY AiRCRAfT DIVISION 
SC.KANLlti.L^ 


VER tin^ni 

9.0 12. so. 00 


ISIltN “ 1 0063S 

IK M ,Nt. 0 I ISICN ■ NODEIII / N U0636 

IH ISlON .Gt. 01 GO 10 ?0 00G37 

00636 

NM.AlIVt MODE NUMBER. IF EVEN* SIGN OF BESSEL FUNCTION WILL BE ♦!. IF 0063V 
01)0, SIGN OF llEiSEL FUNCTION WILL BE -1. 00660 

00661 

IK I M / 2 » • 2 .to. H I ISIGN * 1 00662 

20 C0NS1 « KHUtU ♦ RPRIME 00663 


calculate BESSLL functions of first and sLCOND kind for KNUII I*RPR1HE 0066S 


CALL OESJl const, M, ENJ, TOL, lERl I 
call besyi const, H, EHY, 1Lk2 » 

IHJ ISItN ♦ LHJ 

LMY « ISICN F EHY 

CALCULAlE CHAKALILKISI IC E-FUNCllON 

EMllCU « EMJ ♦ OMUU) • EMV 

60 CUM INUE 
WV9 KLIUKN 
ENU 

fUNcnUN FAL2IP IFUNCT, AL , BR • TOL* ROOT, ITER, YY I 

LUKRLSPliNDS TO DtO VERSION (FALSIE) AHIHJHENT LIST AS FOLLOWS (THIS 
FUR INTERNAL PURPUStS ONLY, IN USE THE TWO ARE INTERCHANGEABLE). 


00666 
0066? 
0066 B 
U066V 
0U6SU 
006M 
006S2 
(I06S3 
006S6 

00655 

00656 
006»i7 
006SH 
006S9 

IS 00660 
00661 


EUNLTIUN falsie (AXR, XXL, XXR, TOL, ROOT, ITER, YYl 


iMlS KOUlINt USLS A CONBINAIIUN OF FALSE POSITION AND BISECTION 
TECmNIUUES to solve FlIK a RUUT (•Hl)Ul*| OF A GIVEN FUNCTION ' 
(•FUNLI*! WHICH HAS ONE ARGUMENT I THE INDl PENDENT VARIABLE I. 

•AL,UR‘ OLEINES THE INTERVAL TO BE SEARCHED. 

The value returned by the FUNLTIUN is FALZIP. FUNCTIFALZIPI ■ ROOT 

THE SEARCH CONTINUES UNTIL TWO SUBSEOUENT GUESSES ARC WllHlN *TOL* 
OE EACH UTHIR, OR UNTIL MUR* ITERATIONS HAVE TAKEN PLACE. 

•YY‘ is RblURNLU AS FUNC 1 (FAL Z I P I , AND SHOULD OE CLOSE TO 'ROai*. 

Ihl iLUiNlUUL HAS ADAPTED ERUM AN ALGO SUBROUTINE APPEARING IN THE 
LUHPUTER journal 12 (19691 — 'EIGENVALUES UF A»X » LAMbl)A»b«X 
WITH band symmetric a and U* by G. PEILRS ♦ J.H. WILKINSON 

LXtlRNAL EUNCI 
real INILKP 
(WRIT ~ 6 

J IS I.OUNI ur ITERATIUNS. 


00662 
00663 
00666 
0066 S 
CC666 
0066? 
0066B 
00669 
006 7U 
00671 
006 72 
00673 
00676 
0U67S 
00676 
006? 7 
006 7 H 
006/9 
0U690 
006H1 
U0nB2 
0U6G3 
006B6 
006(tS 

U06B6 

006(1? 


II 



PKATT t HHJTNty AIRCRAFT DIVISION VER 07/25/77 

SC.PANLIB.14 9,0 12.50.00 


C 

C 


c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 


c 

c 


c 

c 


c 

f. 


c 

c. 


L 

c 


I J »■ 0 
A <• AL 
H <- BH 


fcVALUAIfc FUNCTION AT LEFT (A» AND RlOHT «0| BRACKETS. 

AF = FUNCT |A> 

BF = FUNCT no 

THE FOLLOWING ITHROUGH STATEMENT 3) DETERMINES IF THE FUNCTION 
OPPUSnt SIGN AT THE ENDPOINTS G1VEJ4_ 

ISH » I 

IF (BF > ROUT) 2* 75 t 3 


00660 

00669 

00690 

00691 

00692 

00693 

00694 

00695 
IS OF 00696 

00697 

00698 

00699 


2 ISW e -I 

3 If HAF - ROUT! ♦ ISM) 50, 60, 85 
SIAIEHEN) 5 INCREMENTS THE COUNTER J 


007U0 
00 10 1 
007U2 

FIRST TIME THROUGH GO TO 50. 00/03 

0Q?04 

00705 

00706 
00/07 
00706 

00709 

00710 

00711 

00712 

00713 

00714 

00715 

00716 

00717 

smallest IABS) value OF FUNCTION. 0071 8 

007 19 
00/<!0 

00721 

00722 

00723 
00/24 
00/25 

00726 

00727 
0072B 
00729 
00/30 

00731 

00732 

00733 
00/34 

00735 
TU THL00736 

00/37 

00736 

00739 

00740 


5 J J ♦ 1 

IF LIFT bracket HAS 'SAME' FUNCTION VALUE AS RIGHT, USE BISECTION. 
UTlilKHlSt, SE! UP INTERPOLATLD POINT FOR POSSIBLE USE. 

If (AliSffAF - OFl/BFI - l.E-51 10, 10, 15 
10 IM1LRP = BISECT 
GU 1U 2U 

15 INILRP * <A*BF - B*AF ♦ IB~A)«HUOT> / IBF-AF) 

I 

IF HlTHlM A lULEKANCE OF THE BRACKET B, MOVE THE INTERPOLATED POINT 
ONE lULLKAKLfc AMAY. 

20 IF ( (AtiSflNURP-Bl/ABSIlNT£RP*B) I -2.*TQLl 22,23,23 

22 INII RP B ♦ 1C - B) / AbS IC - B) • TOL 

Sti A-B (B IS ALWAYS THE POINT WITH 

23 A - B 
AF •- BF 


USE PUlNt CLOSEST TU B IIMTERP OR BlSECTl AS NEH B AND EVALUATE BF. 
IF (IINTERP - BISECII ♦ IB - INTERPI) 30, 25, 25 
B » INTeKP 
CO TU j 5 
30 b - BlbfC! 

35 BF = FUNCT IB) 

BFHK = bF - ROOT 

IF CF IS ON THE same SIOE OF THE ROOT AS BF, LET POINT C « POINT A. 
40 IF ULF - ROUT) * BFMR) 55, 75, 50 
50 C - A 
CF - A/ 

IF CF IS CLOSFR (AUS) TO ROOT THAN BF, SWITCH POINTS B AND C. 

IN ANY CASE, B ANU C ARE IHE TWO BRACKETS. ALSU OF IS CLOSER 
ROUT Than CF IS. 

55 IF lAUblBF - KUOTI - AHSICF - ROOT)) 60, 60, 57 
57 A ^ U 
AF ^ OF 


I I 


page is 
POOR QUALITY 


PRATT C MHlTNfcV AIRCRAFT DIVISION VER 07/;$S/T7 

SC.PaNL1B.L4 9.0 12.S0.00 


B ■ C 00741 

BF •> CF 00742 

C « A 00743 

CF « AF 00744 

^ 0074S 

C SfcT UP BISECTION POINT. IF CLOSE ENOUGH* FINISH UP» OTHERWISE 60 00746 

C BACK IF iteration COUNT OOESN*T EXCEED HAXIHUH. 00747 

60 BISECT » (B ♦ Cl / 2. 0074B 

IF { IABS(BlSECT>B)/AB$iaiSECT>6l> -2.4T0LI 75.6S»69 00749 

6S IF CJ - ITfcRI S, 70» 70 00750 

70 WRlTb(lHRIT,1000)J*b.BF,CtCF 00751 

1000 format Uho/// 30X, »1N FALZIP* AFTER' » lA. • ITERATIONS' // 00752 

1 loX, 'BRACKtT I = '. GI5.8, 5x, 'FUNCTION * ', 615.9/ 00753 

2 10X» 'BRACKET 2 * GI5.B» 6X* 'FUNCTION b '* 615.8/ 00754 

3 5X. 'UHACKET 1 HAS RLIURNEO AS RESULT. '| 00755 

75 FALZIP l! 00756 

TY = BF 00757 

KETUKN 0075B 

80 FALZIP B A 00759 

YY B AF 00760 

RETURN 00761 

85 HR11 E(lMRlTtllOO)ROOT«A«AFtBtBF 00762 

1100 FOkHAl J'0***1N FALZIP* ROOT GIVEN (=', 615.8* 'I DION"T FALL BET00763 
iHtL'N VALUES OF FUNCTION AT BRACKETS GIVEN****/ 00764 

2 lOX* 'BRACKET 1 * '* G15.8* 5X* 'FUNCTION c t, gIS.8 / ' 00765 

3 lOA* 'BKACKtT 2 “ G15.8* 5X* 'FUNCTION * •* 615.8 / 00766 

4 50X* 'TERMINATING RUN* I 00767 

STOP 00768 


ENO OOY69 

subroutine INPUTC TEND I 00770 

t 00771 

C THIS SUBRUlillNE INPUTS THE DATA REQUIRED FOR THE EXECUTION OF A CASE 00772 
- ’ 00773 

DIMENSION MXH00EI5OT* ISAVE(50*31 00774 

COMMON /OFAULT/ NLUC* HTR. UR, IMX, FRO* M|50l* HUSI50I* 10IRI50 1*00775 

1 CONDNO* lEHU* XHIN* XMAX* DELX* DELTH* 00776 

2 OELR, NO* XM(50(» KHI50), THHI50I* IPT* NSRH* 00777 

3 NVAL, IXPNT* NPRUB, MPROB, NTAB* NST* SPEED* THIN, 00778 

4 7MAX 00779 

COMMON /CNSTNT/ NHIRES* NHOUES* sigma* B, MX* FREQ* A* OMEGA 007B0 

COMHUN /NUDES/ NUDE 150 1, HU 1 50 1, IHAVEI50I OOIBI 

COMMON /BhSSL/ D0MU5). PI 00782 

COMMON /ANGLES/ DEGRAD* KADDEG ^ 007B3 

common /EMUS/ EMIM 50,501, lEHPRT 00784 

COMMON /BOUNDS/ XuUUNO < 2,50 1 * RBOUND 12 , 50 ) * 1BOUN0(2 * 501 * XA» XB, 00 /65 

1 THMIN* 1HHAX* XLin, RLIH, THLIM 00/66 

COMMON /MIKES/ NOTH* X(50l* RI50), THETA(50| 00/07 

COMMON /search/ NVAR* IPRNT* NSRCH* NEVAL, lEXPNT* NTABLE, NSTAT 00789 
COMMON /PKDBES/ IHALL, I SAME* NPROUE * HNPKOB 0070 V 

COMMON /LUNDTN/ GUUUCN 00/90 

t 00791 

NAMILIST /INDATA/ NUiC, HTH* DR, EMX, FRU, M* MUS, lOlR* CONDNO* 00/92 

I UMU, XMIN* XMAX, tlELX* DELTH* UELR* 00/93 


‘IS 




»»RATT & kHlTNEV AIRCRAFT DIVISION 
SC.PANLIb.L'i 


VER 

9.0 


07/25/77 

12.50.00 


C 

c 

c 


2 N0« XN( RN« THM, IPTi NSRHt NVAL, IXPNT* 

^ NPROB._KEROB^ NTAB, NSTt SPEED, TMIN, THAX 

lENO s 0 

KtA0(5,lN0ATA,EN0&9998l 


CHANGE INPUT UNITS TO INTERNAL UNITS 


20 


AO 

C 

t. UROkK 
c each 
c 


GOGOCN 
NNIKES 
NHUUES 
SIGMA 
b 

MX 
FRtU 
ItMPRT 
XA 
XB 
XLIH 
REIN 
THLIM 
NMM 
IPKNT 
NShCH 
NLVAL 
1 LXPNT 
NPKllBE 
HM'KllB 
NTAULE 
NSTAT 
thHlN 
THMAX 
A 

DU 20 
Hiiok in - 

Muin B 

1 HAVE It) > 

CUNI INUE 
DO AO I^ltNDIM 

xm c 

Kin 

I Ml 1 A I n > 

CGN) INUE 


E CONONO 
S NLOC 
= NUOC 
» HIR 
■ OK 

* EMX 
= FKQ 

* lEMU 
s XM!N 
B XHAX 
B UELX 
« 01 LK 

« OECRAO ♦ OELTH 
B NO 

* IHT 
B NSRH 
B NVAL 
a IXPNT 
B NPKim 
B MPKUB 
B NTAO 
B NST 

B 01 GRAD • TNIN 
B OIGKAO • THAX 
B SPEED 
1b1,NNU0ES 

B Mil) 

B Musin 

* loiRin 


xHi n 

HMl I ) 
THMI I I 


* OEGRAO 


MiiOt NUMBERS FROM LARGEST TO SMALLEST. ALSO ORDER RAOIALS 
Moot NUMBER EROH LARGEST To SMALLEST 


IM NMUOES .EE. 1 I 
lUJ 50 Ul,20 
HXNUliMt) • 1 
50 UNT INUE 

GU 220 1B2, NMUOES 

e MXMODEII) 


GO TO 280 


00 79 A 

00795 

00796 

00797 
0079U 
OOIV9 
00000 
OObOl 
00802 
00803 
OOBOA 

00805 

00806 
00807 
ooeoR 

00809 
00510 
00611 
00812 
00813 
OOblA 
00615 
00516 
, 00817 

00818 

00819 

00820 
00621 
00522 
00823 
00B2A 
0082 5 

' 00826 
00627 
0052 5 
00829 
00B3U 

00831 

00832 
00633 
0053A 

00835 

00836 
008 37 

MlTHlNG0b38 
008 39 
00590 
008A1 
008 A? 
008A3 
006AA 
008A5 
005a6 


n n rt r. r> n 


1 


OK'KJINAl. I'ACK 
Oi? POOK qi;autv 


SC^plNLlU^jr"’' DIVISION 

-0 !m‘ i5?fi 

80 DU 100 J=2ll ‘ 60 TO 120 

rl** » 1 - 4 ♦ 2 

inn ^**800^441 B HXHODt 144-1 » 

100 CONUNUL 

HXHUntm B I 


VER 

«.o 


120 


140 

160 


loo 


OU 200 K«2,l 

c MXHQDEIKI 

IM lAbSl Moot Ml » - UBSI HOOEINEXTI 

hjp W 1 I ^ 


220 


IH MU(1» .LE 
J1 

OU IBO 4=41,1 
JJ 

MXMUtit<J4) 
ttjNT IMJL 
HXmiOL(K) I 


HUINEXTI 
K ♦ 1 


60 To 


1 

200 


1 200, 140, 160 


1 - 4 ♦ 41 
MXH0DU44-1I 


200 

220 


CdNt INiiL 
MXfiOOUl 1 
CUM INut 


CO TO 220 


00 2aO 1=1, 

MAX 

liAVt 11,1) 
liAVt n,2) 
liAVUl,3) 
240 COM INUE 

00 260 1 = 1 , 
MCOl 1 1 ) 

NUI ) > 

IWAVH I I 
260 CUN) INUE 


* 1 


NNUObS 

“ MXHOOEIl) 
MODUHAX) 
« NUlHAX) 

* iWAVtCMAXl 

NHODES 

*= iSAvui.n 
= lSAVEfI,2» 
1SAVE(1,31 


the M«S, Tno nude numbers ARE E«UAt, SET ISAHE INOICATUR TO 1 

9n tfT t c Au ( 


?«0 ISAMl B 0 

IM MUOEII) .Eg. NO0EI2) 1 

CAECULAIL RADIAN EKEQUENCY 

^ 0ML6A B 2. • PI • pRtg 

C 2ER0 tUURDlNATE HOUND AKRAVS 

OU 320 4b1,nlUC 
0(i 300 1 = 1,2 
XIUIUNDCI,4» o 0.0 
Rt'OilNUI 1 ,4) . 0.0 

IHOUMH 1,4) = 0.0 

-100 CCTillNUL 


ISANE 


I 


0T/25/T7 

12.S0.00 


00847 

00848 

00844 

00850 

00851 

00852 

00853 

00854 

00855 

00856 
0085 7 

00858 

00854 

00860 

00861 

00862 

008t>3 

00864 

00865 
00666 

00867 

00868 
00869 
00670 

00871 

00872 
0097 3 

00874 

00875 

00876 

00877 
008)8 

00879 

00880 
00881 
00882 
00883 
UU8B4 
008E5 
00886 
0088 7 
00888 
OOHbV 

00890 

00891 

00892 
00H9 J 

00894 

00895 
0C896 
00897 
00999 
00899 




' I 


»• 


r n 


PKATT L WHITMI^V AIRCRAFT DIVISION 
SC.PANUU.L4 


VER 

«*0 


01/ii/ii 

12.d0.00 


320 CONTINUE 00900 

00 TO 9999 00901 

C 00902 

C END UF DATA SET 00903 

C 0090 A 

9990 lENO • 1 0096S 

9999 RETURN- 00906 

END 00907 

SUuROuTINE KMUCALI VALUEt OEETA, KHUt RlbHT 1 00900 

C 00909 

C THIS SUUKUUTINE CALCULATES THE-CUARAC1 ERISTIC NUMBERi KNU 00910 

00911 

EXTERNAL BESL1» BESL2 00912 

CUHMCiN /CNSTNT/ DUH1I2I. SIGNA. DUM2IS1 00913 

REAL XHUt LEFT 0091A 

30 IPLliS X 0 0091b 

IhlNUS B 0 00916 

36 1F( SlGHA t dOt 90« 60 00917 

C 009 IB 

90 KHU > 0ESL2I VALUE 1 00919 

GO TO 60 00920 

60 KHU B BLSLK value I 00921 

C 00922 

60 1F( KHU 1 80. 66. 70 00923 

66 RIGHT « VALUE ' 00929 

GO TO 130 0092b 

70 IPLUb B I 00926 

GO to 90 00927 

80 IHINUS B 1 0092U 

C 00929 

C OETERHINL IF LEFT AND RIGHT BRACKETS HAVE BEEN FOUND. 00930 

C 00931 

90 IE I IPLUS .EO. 1 .AND. IHINUS .EQ. II GO TO 100 ' 00932 

00933 

BRACKETS NOT FOUND. RECYLL- 00939 

00936 

VALUSV B VALUE 00936 

VALUE B DELIA * VALUE 00937 

GO TO 36 00939 

C 00939 

C BRACKETS FOUND. CALCULATE KHU 00990 

t 009A1 

100 LIFT B VALUSV 00992 

KlbHT B VALUE 00993 

IFI SlGHA I 110. 120. 110 00999 

110 KHU B eaLZIPI BESLlf LEFT. RIGHT. .001* 0.0* 76* VV I 009A6 

GU TO 130 00996 

L20 KHU B FALZIPI 8ESL2t LEFT. RIGHT. .001* 0.0* 75. VY I 00997 

130 KEIURN 0099B 

I NO 00999 

subroutine KUCAL 00960 

L 00961 

(. THIS SUBROUTINE LALCUEATES THE CHARACTERISTIC NUMBERS KHU AND DNU 00962 


ORIGINAL PAGE IS 
OP POOR QUAUXB 




07/2i/TT 

12.&0.00 


OWICW.Al. I- , 
POOK 


9.0 


COMMON /Moots/ M0Dfc|50l» MUISOI. 1MAVEI50I 
COHMLjM /KOHU/ KMUI50). QMU<SO) 

COMMON /btSSL/ ISlGNi JStGN. OELKNU. TOL H pi 

00 ioo i=i,NMunes 

CALCULATE UKOtK FOR BESSEL FUNCTION EVALUAJ40N 
lA&St HOOEIU ) 


M 

IFI N 
ISIGN 
JSIGN 
BKAKTL 


>NE. 0 ) 

“ I 
= -I 

« .1 


C 

c 

c 

c 


10 ISIGN 
JSIGN 
ERAKTL 
IFI ISIGN 


* MOOEdl / H 
ISIGN 
•• M 

•G£. 0 I 


GO TO 10 


GO TO 20 


CO TO 20 




c 

c 

c 

c 

c 


IFI I M 
20 NtIMMUS 


/ 2 » • 2 .to. M I 
- Mum ♦ 1 


ISIGN 


C 

c 

c 

c 


•Hi «SUN m SINUU»»“™s 'onJ MU.flOH OtflNINC 

KMUPkM . 0.0 ' 

00 ^G J:>l(NUHMiJS 

IFI M .LU. 0 .And. J .to. 1 l go to 

OELKHU, KMUPRM, BRakTR ) 

BKAKTL s BKAKTR 

*0 CONTINUE 

FMUIll « KMUHKM 


GO* 60t 80 


tFI SIGMA » 

60 «MUm * 0.0 

BO IM KMUIII » 

NO C«L »ISJ, K„,p.N. N-JSIU., IMNl.^Ilt ll”, 

'S>GN ♦ JSIGN ♦ EMMl 
CALL BESJI Kmiimkm, H, £MJ, TOL, IER/ ) 

CAll. htsn MliJPHM* M. EMY, ILRJ ) 

•- ISIGN • tMj 


00953 

00954 

00955 

00956 

00957 

00958 

00959 

00960 

00961 

00962 

00963 

00964 

00965 

00966 

00967 

00968 

00969 

00970 

00971 

00972 

00973 

00974 

00975 

00976 

00977 

00978 

00979 

00980 

00981 
100962 

00983 

00984 
0096 5' 
00966 
0098 7 
00988 
00969 

00990 

00991 

00992 

00993 
>00994 

00995 

00996 

00997 

00998 

00999 
GIOOO 
OIUUI 
01002 

01003 

01004 

01005 


•><> 


tv 


I 


PHATT i. WHIINI Y 
it.PANl 


AlRtRAPT OlViSlUN 


VtR 

9.0 


ulsyi 


I MV 
( ALL 
L VMl 
A 

h 

Sli.MAK 
LALL liLSJI 


HtSj< 

PtSYl 

IUiY< 


100 

•#999 


I ALL 
LAt L 
CALL 
LMMl 
I YMl 
I MJ 
1 MY 
L 
U 

wHlli i I 
I i.Nl IMOL 
Kl IIIKN 
1 Nil 

SUPKI-Ltl INI 


k 1SU.N • LMY 
KHUPKM, M-JSlf.N. 

U IMl.N * JSU*N 

LMMl - 1 M • JSION • tHJ 
.. IVMI - < M • 

= SU.MA ♦ KMlIt’RH 
S.H.MAK, M-JSUN. 

Sll.HAK, M, IMJ, mit UR6 J 
Sll.MAA* Mt • 

MI.MAK, M JSH.N. 

^ ISUN • JiU.N • IHMI 
.. iSUN * JMtJJ • tYMl 
= isn-N ♦ IMJ 

t tSll.N ♦ IMV . ^ I 

t IP.MI - I M * JSlON ♦ iMJ 

. lYMl - t M • JSK.M * 


kmuprm 

KMUPRM 


UR6 » 


^ - 1 A ♦ C ♦ 


Sli.MA » / « 0 ♦ « 


SIGMAK 
SIGMAK 
* SIGMA 


print 


c THIS SimROOllNl prims input 


and CALCULAim VAlUfcS 


PRLO» A* OMtGA 


I MU I *>0fS0l t It Ml Kl 

NUlH, M*> 0 ». R<^ 0 >« 
NVAK, ll’KNIt NSKCU. 

ilMMli'. /WAVlNti/ KAf'lll 

/KyMO/ KHUlt.O), UHOl 01 

him MX, RMO 
Cli:SPL» X KX • KXANGl 


{. I MMl N 
CiiPHliN 

1 

( IIMMIN 
tOMMliN 
I N 


/ANULI S/ 
/UlUINUi/ 

/IMUS/ 
/hlM S/ 
/SI ARLH/ 


1IU lAlSOI 
NkWAlt ILXPNli 


XA. XB 


NTABLLt N^1*.T 


L 

t 

C 


UiNVl Kl 


INllKNAl llNilS in nl'lPOT ONUS 


1 HI I H 
lliMlN 
IMMAX 
LU1 “lO 
liU .'0 


i‘0 


AO 


1 KlIilNlU I f*l 
I I NI INOl 
1 1(1 1 AlJ > 

I I NI INlIt 


>- RALHUG 
. KAlHil G 
^ KAIHH G 

J^UNMlMS 

l-l..’ 

RAUOLG 


- HAODI G ♦ 


IHLlH 

IHMIN 

lllMAX 


1tUillN0ll«J> 
Ihl t At Jt 


10 


UlNVIHl ANY NtlAUVL AN(.U S 

IMI Ahl.POSt IIIHA. NMIKLS 


PtlSlUVt ANCLtS POR PRIN1ING 


0I/2^/T7 

12.^0•00 


0100 b 
01007 
OlOOB 
OlOOY 
OlOlO 

olon 
0)012 
oloia 
D101<> 
OlOlb 
0101b 
0101 1 
OLOIH 
0101*7 
01020 
0102 1 
01022 
01023 
0102 A 
0102 b 
0102b 
tilUi’T 
01O2B 
(jl029 

01030 

01031 

010i2 

01033 
0103A 
0103b 
01036 
010 3/ 
OLOjB 
oIoav 
OlOAO 
OlOAl 
010A2 
010 A 3 
0 ioaa 
OlOAb 
OlOAb 
010 a 7 
OIOAB 
010A9 
OlU'iU 

C 1 OS 1 
01oS2 
OluS3 
OlOSA 
OlOSS 
0 lo'ib 

Olob f 

ulusn 


■A\ 


I 




PRATT & KHITNEV AIRCRAFT DIVISION 
SC.PANLIb.t.4 


VER 

9.0 


C PRINT INPUT 
C 

WRnb{6(9000l 

4000 format I IHlt T44| MICROPHONE LOCATION COMPUTER PROGRAM 

1) 

MR11b|6, 90131 

9013 FUKMATI //, T66i *... INPUT VARIABLES I 

HHllblb^VOOn NHlKbS* NMCjObS 

9001 FURMAH //. TS, ‘NUMBER OF MICROPHONE LOCATIONS * 'f I2i T5li 

I'MJrtBtR OF (HOOE.MUl SETS » 12 I 

NKlTtl6.y002) 

9002 FOKMAH //. IX, •... INPUT MOOES .»»•» //• T5» 
1‘LJKCUMFtKbNTlAL*, T34, »RA01AL*, TAT, *MAVE», 


2MblRS lev, 'KMU*, /, Tl&f 'MODE NUMBER*, T34, 
3*1M<1CATUK*, 162, ‘REAL*, T71, •IMAGINARY*, f ) 

DO BO t-l,NMO0ES 

KXANGL = RAODEG * KX ( I) 

HKITE 16,9003) 1, KODEU). MU<U, IMAVEII), KXANGL, 
9003 FUFHATT 6X, 12, IIX, 14, 13X, 12, UiX, 12, 9X, FIO< 
ISX, Flu. 9 ) 

BO CONTINUE 

Wknbl6,9005) SIGMA, B, MX, FREQ, A, OMEGA 


•MODE*, T14, 

T61, *AX1AL HAVE 
•ORDER*, T47, 


KMUll) 
4, IX, 


FLO. 4, 


9003 FORMAT { 

1 • HUB 
2T0‘*, 

3T42, 'SPEbO 


///, IX, *... TEST geometry AND CONDITIONS ...*, 
/ UP RATIO = *, F8.3, T42, ‘OUTER RADIUS OF DUCT - 
•AXIAL MACM NUMBER = *, FB.3, /, T&, 'FREQUENCY * • 

F9.2, TB4, 'RADIAN FREQUENCY • 


IJF St UNO *, 


//, T5, 

*• F9.3. 

) FI0.3, 
>, F12.3 


//, T5, 


T20, 
•MAXIMUM R*, 


HKirU6,9006) NVAR, NSKCH, NEVAL 

9006 FO'.MAH ///, IX, SIQLMASTIC SEARCH V*LUES . 

l*f.UMbhk LF VARIABLES = *, 12, T42, 'MAXIMUM NUMBER UF SEARCHES 
2, 13, 1B4, 'NUMbbK UF EVALUATIONS PEk SEARCH » *, 15 1 

RK11L(u,9U07> 

9007 FOKMAH //, 16, *M ICROPHUNb • , T69, 'PARAMETER BOUNDS',/, 
I'MINIMIJH X', VX, 'MAXIMUM X', 6X, 'MINIMUM R*, 4X, 

26X, 'MINIMUM ThEIA', 4X, 'MAXIMUM THETA*, / I 

DU IDU J^2,NMlK.ES 

HkIIL(6,9ua8) J, I XBUUNDt 1,JT ,1=1,2 I, \ RBOUNDU ,3 1 , 1* 1 , 2 I, 
1 I TBCJUNDU.JI, 1 = 1*2 ) 

900B FOKMAH 6X* 12, 3X, 2 HX ,F 7, 3,6X,F7. 3) , 9X, F7.3, lOX , FT. 3 I 
100 CONTINUE 

L 

C PKINT HILKOPHUNE LOCATIONS 


HHlTHb, 90001 
HklTM6,9009> 

9009 FORMAT t //, IX, 

I'MlChUPhUNt', T25, 

DO 120 Ul.NHlKES 

WkllL(6,9UlU» 1, XU), Rin, THETAU) 

9010 MA<HAT( 4X, 12, 9X, 3(F12.6,6X) ) 

120 CtNT INUE 


, CALCULATED MICROPHONE LOCATIONS 
»X», T42, »R', T57, 'THETA*, / » 


// 


PRINT LHAHACTEhlSTlC E-FUNCTlUN VALUES IF REQUESTED 


07/25/77 

12.50.00 

01059 

01060 
01061 

***•01062 

01063 

01064 
01066 
01066 
01067 

01065 

01069 

01070 
NU01071 

01072 

01073 
0107V 

01075 

01076 

01077 
0)078 
01079 
01050 
01081 
010B2 
01063 

lOlObV 
01086 
01056 
» *01087 
oloaa 

01089 

01090 

01091 

01092 

01093 

01094 

01095 

01096 

01097 
01095 

01099 

01100 
01101 
01102 

T5, 01103 

01104 

01105 

01106 
01107 

oiloa 

01109 
OHIO 
01111 






1‘RATT t WMITNEV aircraft OlViSIllN 
SC.PANlIR.tA 


VtR 

9.0 


tM UHPHi *lt. 
WRItt(6*90UJ I 


0 ) 

Itl^ltNNODES 


GO TO 9999 


» 


IXi 


VOll HIKMAll ////t 
UllMIN VACUl S) 

1^)« f I 

lUJ 1-0 J^^UNMIKCS 

hKHUht'^OWJ J» * fcHUU f Jl»l = l»NMOl)fcS 
KIKHAU AXt !?• lS(lXiF7.3) I 

COM INUl 


CMARACTtRlSnC E-fUNtllON VALUES lUMAL 
//, IXt * 100 X 1 ION* • 168 » *H 00 tS*i /i »*l 


9012 

lAO 


9999 KLUIKN 

tMO ^ 

StiORUuTlNfc SlUCHt FCT* 1 N 0 | WORK ) 

UmtNSlllN INDUI* HORKm 

UilMuI^/lNMt/ NPARM. NSfcRCH. NEVAU. '♦'ABLE, NSTAT, ISKO, lOlST, 
ItXPNT, IPRINT, ITYPh, LT, JT, ICNVRG, TOL 
/INOEX/ DUM 112 ), INOEXA, INDEXS, 1 NDEX 6 , INDEX? 

/SfcLO/ JStED 


i 


LUNHUN 


C 

c 

c 

c 

c 

c 


INlTlALl/l PARAMEILRS tOK SEARCH 
CALL INllALC WORK, IND » 


SEAKUl UlOP 

UU 120 
JSi to 


l*^l,NStRCH 

ISEED 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


GENIRAU STAllSlICS FUR SEARCH I 
CALL SI AH WORK, 1 > 
tVAHlAlilCN LOUP 

DO 100 J*^lfNtVAL 

GENERA1E RANDOH PARAHETER VALUES 
CALL RANtUlHl WORK I 

IVALUAIL M1NC110N MlTH RANOOH PARAHETER VALUES 
3 FCTI 0* NURKI INDEX7«1 1, lb 


UP 

IP HINCllUN VALUE 
NEXl i VACUA! IIIN 


10P> IS HEANINGEUL UPDATE G TABLE. IF NOT GO O'N TO 


loo 


1E( IE .01. 0 > 
LAIL UPKAU. ( KURK, 
CUNT INUl 


GO TU 100 


I.P 1 


PHlNl Kt SULIS UF SEARCH I IF REUUFSTtU 


07/'26/77 
12. SO. 00 

01112 
01113 
SOLO 11 19 
01119 
01116 
01117 
OHIO 

01119 

01120 
01121 
0U22 
01123 
01129 
0112 '> 
01126 
01127 
0112U 

01129 

01130 
on3i 

01132 

01133 
CI139 
01139 

01136 

01137 

OlUa 

01139 

01190 

01191 

01192 

01193 
01199 
01199 
0119b 
0ll9l 
0119B 
01199 

01190 

01191 

01192 

01193 
01199 
01199 
0119b 
01197 
0119U 
01199 
OllhO 
Ollbl 
0Hb2 
01163 
01169 


OiViCVM'Al- r li’ 

jOt: ('t A: ' ' ■ 


PRAll L WHITNEY AIRCRAFT DIVISION 
SC.PANUlb.L4 


IPR1N1 


CALL PR1NTH WORK, 1 


U ST FUR CUNVLRGtNCE. IF CONVERGtO, RETURN WITH VALUES OF VARIATES 
IN The first NPARM LLEMENTS of work array, if not I CO ON TO NEXT 

SEAKU). 

CALL CUNVRGt WORK 1 

IM ICNVRG .G1 . 0 » GO 10 140 

KiitALL ® fCTl 1, WORMINDEX?*!!, U 1 

120 CljNT INUL 


;1NT results OF LAST SEARCH IF REQUESTED 


1-1 
. 1 > 


lAO 1F< IPRlNT .EQ. 1 > CA 

IbtST = 1F1K< W0RMIN0EX4*U 

DU J60 1*1, NPARM 

jhDEX •= 1N0EX6 ♦ NTABLE * I 

WORK ( 1 ) HURKC JNDEX* IBEST ) 

16 U LUNTINUL 

1NUI61 = ISEEO 

9VY9 return 
END 

SlIFKtJUlINE INITALC WORK, IND 1 


CALL PRINTTI WORK, 

II ♦ .3 I 


^UBKIIUT 


INL initializes PARAMETERS FOR THE STOCHASTIC SEARCH 


DIMENSIUN HOHRUI, INDIll , 

tD^v^UN /INITL/ NPARM, NSERCH, NEVAL, NTABLE, NSTAT, ISEEO, IDIST 
I lEXPNT, IPRlNT, ITVPE, LT, JT, ICNVRG, TOL 

LllMHUN /INDEX/ IN0LK2, 1N0EX3, 1N0EX4, IN0EX5, 1NDEX6, INDEX? 
CUMMIN /LOOPS/ NLOOP 

NPARM = INDllI 

NStRCH “ 1NDI2I 

NEVAL = INDOI 

NIAULE “ 1NU14) 

NSTAT - INlMSl 

ISLEU » IN0I6) 

101 ST » IND(I) 

IIXPNT = INDIHI 

IFRINT * 1NU191 

HYPE * INDllUl 

INDEX? * ? * NPARM 

INDt X3 ». 3 ♦ NPARM 

lfll'EX4 * 4 * NPARM 

Jfi.U xtJ = INDl X4 ♦ NTABLE 

iN|(i.X6 = INDtXb ♦ NTABLE 

INlHX? - 1NUIX6 ♦ NTABLE * NPARM 

NLUUP * NPARM / ? 

IM T * NLUOP .Nl . NPARM I NLOOP ■ 1 ♦ NIOOP 


NLOOP 


07/2S/T7 

12.30.00 

01163 

01166 

0116? 

01I6B 

0116V 

01170 

oim 

01172 

01173 
OH/4 

01175 

01176 
0117 7 
0117B 
01179 
OLIBO 
OllBl 
011B2 
01183 
D11B4 
01163 
01186 
01167 
OllBB 
OllbV 
OllVO 
OllVl 
011V2 
011V3 
011V4 

T, 0UV3 
0UV6 
01197 
01196 
01199 
01?U0 
01201 
0120 ? 

01203 

01204 
0J2US 
01206 
01207 
0120H 

01209 

01210 
01211 
01212 

01213 

01214 
01213 
Ul2l6 
01217 


r 


oof' no,^r'n or. n 


I 


QUAXJl'Sj 


PKA1T L HHIINtY AlRCKAfT DIVISION 
Stit'ANtlb.L'i 


WER 07/25/77 

9.0 12.50.00 


LT 


nTable 



JT 

& 

NTABLE 



1 CNVKO s 

0 



TUL 

ft 

10. ** 

FLOAT 1 

lEXPNT 1 

in 

NTABLE .ll 

. 1 ) 


NTABLE > 2 

in 

NSTAT .LE. 

1 1 


NSTAT « 2 

1F( 

NSTAT ,&T. 

NTable 

> 

NSTAT *■ NTABLE 

IF( 

ISLU) .LE. 

0 1 


1st ED a 1 


C SET K, G» AND (> TAULES 
C 

00 20 Isl.NTAhEE 
WUKM1N0LX4+1) » 1 
20 CONI INU! 

00 AO l-l»NTAtiLE 
WOkAUNoLXbtl) 0.0 

AO CONI INUt 

LtMLNl * NIABLE • NPARM 

OU 60 l^liLEMbNT 
WOkK UNlilXo^l t s 0.0 
60 CUN1 INUL 
VV99 KbIURN 

1 Nu 

SubiUiUTINE S1ATI HURKi ISERCH t 

This suBRooiiNt genekates the statistics required for The search 

UlMbNSlUN MDRA(I) 

CUhMoN /INITL/ NPARH, DUH1UI» NIABLE, NSTAT , 0UH2I9I 
LCMHliN /INUtX/ lN0bX2, 1N0EX3, INOEXA, 1NUEX5, 1N0EX6, INOEXT 

IM ISERCH .Cl. I » CO TO 40 ’ 

FIRST SEARCH - CALCULATE MEAN AND STANDARD DEVIATION FROH INPUT 
PARAMEILK BOUNDS 

OU 20 l-l,NPARH 
BL - HOKKin 

6R - MURK(NPARH«n 

HURK(INDLX24n a .5 * ( BR 4 BL I 
Wt.RK(INUbX3*U a .5 * I BR - BL I 
20 CONTINUE 

GO TO 9999 

OEIERMINE MEAN ANO STANDARD DEVIATION VALUES FRQN (* TABLE 
40 DO 160 Jal, NPARM 

JNDEX “ 1N0EX6 ♦ NT ABLE • 1 J - I I 

PAVG *= HURKI JNDEX41 1 

OU 120 Ia2,NSTAl 

PAVC “ HURXlJNDEX*n •* PAVG 

120 CONI INUL 


0121U 

01219 

01220 
01221 
01222 

01223 

01224 

01225 

01226 
01227 
0122H 
01229 
C123li 

01231 

01232 

01233 

01234 

01235 

01236 

01237 
01239 
01239 
0I24U 

01241 

01242 

01243 
01294 

01245 

01246 
Cli:47 
0124B 
01299 

01250 

01251 

01252 

01253 

01254 
01235 

01256 

01257 
0125H 

01259 

01260 
01261 
01262 

01263 

01264 

01265 

01266 
01267 
01265 

01269 

01270 


r> n n mo 


- rooii 


' ^ ir. 

j'l\' 


PRaTT & WHITNEY AIRCRAFT DIVISION 
SC.PANUD.I^ 


VER 07/25/7T 

9.0 12. SO. 00 


PAVU B PAVC / NSTAT 01271 

11 •= IFIXI M0RM1N0EX4«J| ♦ .3 I 01272 

WCjRK(1N0EX2«J) s waRKiJNUEX^ll) 01273 

PSIGHA s 0.0 01274 

UQ 140 UltNSTAT 0l27d 

PSlOHA B ( waRKUNObX^ll > PAVO ) *9 2 * PSIGNA 012/6 

140 CONllHUE 01277 

HUI«K(lNDEX3«J) a $URT( PSl&NA / | NSTAT - I I 1 01278 

1'60 CUNUNUl 01279 

C 01280 

9999 RLTURN 01281 

tNU 01282 

SUURUU11NE RANDOHt WORK ) 012&3 

C 0)284 

C THIS SUbKUOTINt GbNbRATbS RANOQN VALUES USED IN THIS EVALUATION 01286 

G 01286 

OlHENSiON WUKKll) 01287 

CUHHUN /INITL/ NPAHN. DUN1I2). NTABLEi NSTATt ISEEO. lOlST. 01288 

1 miH2i7) 01289 

COHMIJN /INUEA/ lNUtX2i INDEXBt 0UH3I3)« 1N0EX7 01290 

COMMON /PHUBbS/ IWALL. I SAME . DUH4(2) 01291 

COMMON /LNSlNT/ 0UH6I3), 8. IIUM6(6» 01292 

COMMON /MIKES/ NOIM, XISO), K160), THETAI501 01293 

COMMON /BOUNUS/ UUH7I301). X8f OUHH(S) , 01294 

C 01296 

DU 340 J^ltNPARH 01296 

1F( ISAME .OT. 01 GO TO 120 01297 

01298 

MODE NUHbER IS NOT THE SAME AS THE FIRST MODE NUMBER 0129V 

01300 

IF( IWALL ) 40( 20 * 60 01301 

20 1F< J .LE. 21 GO TO 280 01302 

W0HK(lNDEX7tJ) a 8 * 01303' 

GO TO 9999 01304 

40 lEl 4 .GT. 1 ) GO TO 280 01306 

H0KKUNDEX7«J) a XB 01306 

GD TO 340 01307 

60 IFI J - 2 ) BO* lOOt 280 01308 

80 H0RK(1NUEX74J) a X(NDIN-1| 01309 

GO TO 340 01310 

100 H0RK(INUEXT+J| a THE7A<N01M-1I 01311 

GO TO 340 01312 

01313 

MODE NUMbER IS THE SAME AS THE FIRST MUUE NUMBER 01314 

0131S 

120 IFf IWALL ) 140, 160, 220 01316 

140 IM J .GT. t ) GO TO 2B0 01317 

H0KKnN0LX7«Jl a XB 01318 

CO TO 340 01319 

160 1F< J - 2 I 280, 180, 200 01320 

180 WCiKK(lNllEX7«J) a 0.0 01:i21 

GO TO 340 01322 

200 WQKKUNDEXTtJ) a 8 01323 


1 -, 


55 


,.raT1 (. WMITNIV AIKCRAM DIVISIUN 

St.H ANl A 


VER 

9.0 


4-.(0 IH .1 - 2 » 

241) WllHK ( INUt A M Jl 

2<»0 HIKK UNIU.X7* J> 


X«N 0 lM-n 

THtlA<N01H-ll 


(iO TO 9999 
240. Z60, 260 

GO TO 3^0 
Gl) TO ^0 


t GlNLKAH HANOUH VARIAOLb 

„0 LJLI. 51V.«< . »U«KlINDtK3.J. 

‘c CH,C^ .».! 1...^ V«U, 0^ « IS KMU.R13 .OURRS. 

c KU.tNlKATL A NIM VAl Ul Ob XX 

^ If I XX - WORKUI I 

30U IM XX - WUKK(NI*AKM«J» > 320, 

y/u wiiKKunm xMJ» * XX 
J<,U t.l'Nl JNUl 
L 

SV'W Rf UIKN 

SOoHliUllNb UI'OAItl WORK, 6P ) 

t m.»L UH..AUS IHt 0 «»,. .« CU«RA». V.LU, Of Of 


lASTl. 

I.UST 


rSm^'^mlu^'llJAKR. OURUfl, OIAOLE. 00K(0>. HVK. •>’ ' 

/.NOL.f 

= 1HX( miKMlNDt XA.NlABLbJ ♦ .3 » , 

= HURMlNntX!fbASl(.I ' OP H0RKUN0EX5*LASTGI 

IH IIYPL .01. 0 > 

^ eilHPAKL OP HITH 1Mb MNIAIUU bUHtNl Of 1Mb G lAlUb 

L , , 9999, 9999, 20 

Ifl Olbbl » 

20 K1 “ LAStG 

HUKKllNI)bX?>*LASlG> ■= GR 

,,J, «. - X . . X.S10 

iuKKljWlXKl • HORRllUOtXXfll 

40 UlNllNlIb 

L OfllKMlNl fUSniON Of OP IM G TABLE 

|i?xn. ^ MimMlNDlX<,*Jl) ♦ .3 I 

. HIIRM INl»tX?»*NtXlOI - tp 

Ml IlVPb .01. 0 » 
ui oil SI 1 
HU IM J1 «Ol. I 1 


GTtST - «» - MORM1N01X5 ,NEXTo» 
100, lou, 60 
Gl> 111 GO 


07/25/TT 

12.40.00 


OI32*t 

0132b 

0l32h 

01327 

01320 

01329 

01330 

01331 

01332 

01333 
0133A 
0133b 
01336 
0133/ 
0133B 
01339 
013AU 
013A1 

01392 

01393 
01349 
0139b 

01396 

01397 
0139 H 
01399 
Ol3bO 
01 3->l 
013b2 
0l3b3 
0l3b9 
013bb 
01396 
013b7 
Ol3bB 
Ol3bV 

01360 

01361 

01362 

01363 
01369 
0136b 
U1366 
01367 
01 3611 

01369 

01370 

01371 

01372 
01313 
01379 
0137b 
01316 




P'nc r>pr> ooo o 


/ POOR QUAIOT 


PKAIT U WHITNEY AIRCRAFT DIVISION 
SC.PANLlb.L4 

JT B 0 

100 IFJ JT .CO. LT-I I 


VER 

9.0 


GO TO 120 
GO TO 160 


C update K (INDEX) ARRAY TO REFLECT ENTRY OF GP INTO THE G TABLE 
C 

120 W0KK(1NDEX4«LTI » W0RRUNDEX4«LT>U 
1F( JT .GE. L(>2 ) 

LT = LT - I 


GO TO 140 
CO TO 120 


140 H0RK(1NDEX4«LT>1) s KT 
160 LT ■ NTABLE 

JT • NIABLE 

V99V RETURN 
END 

SUbKOUTLNE CONVRGI WORK ) 

THIS SUBROUTINE TESTS FOR THE CONVERGENCE OF THE INDEPENDENT PARA- 
H t T L R S 

DlHtNSUm WOKKIil 

CUMMUN /INITL/ NPARNt DUMK2), NTABLE* DUH2(6)* iCNVRG* TOL 
CGMMUN /INDEX/ DUM3(4), 1NDEX6* INDEX? 

COMMON /PROBES/ IwALL* ISAME. 0UM4(2) 


CALCULATE MINIHUH AND MAXIMUK P VALUES 

3S JNDEX B 1NDEX6 * NTABLE • I I - I ) 

PMIN B WURKI JNDEX41) 

PHAX B PHIN 

DO BO JB2, NIABLE 
P B WOKKfJNDEXi^J) 

IF( P ~ PMIN I 40, 45, 45 


07/2S/7? 

12.50.00 

01377 

0137B 

01379 

01380 

01381 

01382 

01383 
01364 
01385 
01385 
01357 
013HH 

01389 

01390 

01391 

01392 

01393 

01394 

01395 
01395 

01397 

01398 

01399 
U14UU 



DD 

U ( 

ICO iBi.NPARM 
ISAME .01. 0 1 

GO TO 15 

01401 

01402 

01403 

MODE 

NUMBER IS NOT EQUAL 

TO THE FIRST MODE NUMBER 

01404 

01405 


IF( 

IWALL ) 

35* 5* 10 

01405 

01407 

5 

1F( 

I .Lb. 2 ) 

GO TO 35 

01408 

to 

1F( 

1 .Lb. 2 ) 

GO TO 100 
GO TO 100 

01409 

01410 

MODE 

NUHbER IS EQUAL TO 

GO TO 35 

FIRST MODE NUMBER 

01411 

01412 

01413 

15 

IFI 

IWALL ) 

20, 25, 30 

01414 

01415 

20 

1F( 

1 .bO. 2 1 

GO TO 100 

01415 

25 

1F( 

1 .Gb. 2 ) 

GO TO 35 
GO lU 100 

Ol4W 

01415 

30 

1F< 

1 .Lb. 2 ) 

GO TO 35 
GO TO 100 

01419 

C142U 


01421 

01422 

01423 

01424 

01425 
01425 
01427 
0142b 
01429 




r>onr% 


^0 PMIN B p 

IF( P - PHAX ) 

60 PMaX a f, 

•»y OilNTlNUt 


VEH 

9.0 


GQ n> eo 
GO* ao* 60 


IM TLiT - TOL ~ » / AHAXil I., aBSI PMAX ♦ PHlN 

100 continue lOOf 100* 9999 

“"du,“JS‘ "" «« «««««.. SU IC«,„KE 

ICNVRO a t 

»999 KMURN 
ENll 

iUt.KUillNt SIVAKt J. C. X, Ifi » 

►«»- ^>t"« . «.»«c o« 


^ lf< J .GT. 0 > 

C NOKHAL 01!klKI|tUTlUN 
K a 0 

20 K a 1 H 

AO CALL RaNO« Ro, la I 

•^1 e RO 

60 CALL HANOI R/, la i 

JH R2 - R1 , 

GO CALL HANOI RL, la 1 
IM R1 - R2 J 
100 CALL RANOI RJ, la ( 

a X 

in H2 - ) 

C • I PK ♦ RO » 

“I- • < l-K ♦ RO > 

C 

C hXPONLNTlAL OISTRIHUTIUN 

160 K ^ Q 

IbO K 

?00 call HANOI K(| 

HI 

220 CALL Hanoi r2, 

•f < K2 - HI ) 

260 LAI.L RANOI HI, 

IH K1 - R2 , 


GO TO 160 

GO TO AO 

BO* 100* loo 
60* 20* 20 

120* lAO* lACL 
GO TO 9999 
GO TO 9999 


1 

♦ K 

GO TO 200 


IH 

1 



HO 




JB 

1 



111 

1 

«A0* 260* 

260 



220* laO* 

lao 


07/25/77 

12.60.00 

01930 

01931 
0X932 
01633 
01639 
01635 

Hul936 

01a37 

Ul<,38 

0l63v 

01660 

01961 

01662 

01693 

01999 

01995 

01696 

U1697 

AN0199R 

01669 

01650 

01961 

0195? 

0165^ 

01656 

01956 

01956 

Ul65V 

01968 

01969 
Ol9bU 
01661 
0166? 
01663 
01666 

01665 

01666 

01967 

01968 
0K69 

01970 

01971 
0197? 
01671 
016/9 
01975 
019 76 
01677 
0 16/11 
016/9 
01980 
01681 
01682 


01 -' • ' 


PRATT C WHIINIV AlRtRAfT DIVISION 
St. pant 1W.L4 


/60 TK 
X 


CALL KANDI K2, lU ) 
IM K? - t » 
l.AiL KANO I R^, IB » 
TK K»’ - .S » 
iUO X • 


K 

fK ♦ KO 

tXPT X-.S*X*X-.S» 


3.^0 X 


C ♦ X 

-< • X 


?H0, ?eO, 160 
300t 300. 320 
bU TO 9999 


9999 KITURN 
(NO 

SUf.RtlUTlNt RANOI VTL, IX ) 


- IX * 6AS39 
IM IV .LI. 0 t 
V»L ^ JY 

“ .A6S6613L-9 ♦ VfL 
IX IV 

9999 hLTlIKN 
(NO 

MmRtXlTlNfc PRINT I I MllRKt IStRCH ) 


IV ■ 21AS4(t3697 ♦ 1 ♦ IV 


07/?S/77 
I2.su. 00 

0IA93 

01AH9 

OlAHS 

0UH6 

UlAb? 

OlAHB 

019119 

OI9Vl> 

01991 

0199.’ 

01993 

01999 
0U9S 
01996 
01‘.97 
019VB 
01999 
UlSUO 
OlSOl 
01902 
01 SO 3 
01S09 
01 SOS 

01506 

01507 
OISOH 
0lS(i9 


C 'sLAVmr""''^ ASSOCIaTLD PARANtltR VALUfcS M.Roisio 

L OISII 

01H.NSK1N WURMIT, rTVPL< 2». KniST(3,2l. 1VARI3I. SDLVUI 

I *sitn. lOisi. oisia 


‘ KXPNT. IPKINI. IIVPI, 0UHH31, ilil ’ 

/'.‘s:;.;; ‘Kp* 

a.MRlN /Mil,/ JMIl) 

UiM.Hl.N /ANOM S/ III OKAO, KAOOl {, 
t'AlA RTYPI / hMHINI. MIHAX 1 / 

I'^Ia T • AmXPd, 9HNINT. ATllAL / 

^ lAIA IVAk / sH X , AlllMlIt 9(1 R / 

^ >n IMKIM .Nt. 1 .im. IPKINT .Nl. 7 I GO TU 20 

C PRINT INPUT PARANlIlRS 
L 

KRmiN.lOUOl MVPI IIIVPIAII. 4 M)1SI4I.|1>IS1«1 |.i-I. 3 I, tOL. 

IOOo\|1kMAT( // 1 NSIHCII, NfcVAi» 4SI Hi 

1000 H1KMAT( //» »»'•. *A *, A9, •HUM IS BilNO StUlOMT*, //, 1.36, 


0 r> I s 
01S16 

oisi r 

OISIB 

01M9 

OIS/O 

UIS2I 

01S2? 

Olf./,3 

0TS29 

OIS/S 

OIS/6 

01S77 

01 !i?B 
01629 
OIS TO 
1) 1 S U 
01932 


A A. A I.^O* 

> .f,,, * •OlSlKIMlHuN IS 01 INI. USIll*, //, T36, OlSll 

' MlUiKANll nm UINVlKl.iNCl *, Tils, •;. •, lilo.9. //. t 36 ills** 

./y t,.. r r ‘ ‘HHIMS in lUIIKHlNANI TABU*, IBs, •, 16. IHS 3 s 

>//, 136, NiKiniR 111 I 1 1 hi NT S USI (. lUMNlRAtl STAHSHlS 16,//01S3S 


r.ih, 


VtR 

».o 


07/2S/7 

l?.^O.Ol 


c 

L 

L 

( 


■' *o.'»o 

HUM KANOllM IMUMWlRt, Ik^, ^tARCM*, TB01S<7 

» t I 

Nl Vt AlyiiAt.<<k . 


SI ANDAM) (UV 1 AI 1 iin< 

- ::::‘r “ - - »*- 

I'l' * s 1-1 ,NRakM 'U 30 

wiNn*m’t ‘ “I’KMiNi.txj*!) 

' IIAOIUG • SDtVf^l 

AM* Tl u 4. A . ... 


PH INI 


lAl'tt ANl, **AKAMtT(K VAlUlS 


7 ) ) 


-oi* 

, *MSuns u>- „a AFt.« . 

►'Hllf (6,10„. , , IVAKIII • . . StARCH* 

;•" - ......A,,:. • "■ ^ , 

AMU ; ♦ .3 I 

H khah sx, , 

‘•t- iiV.JiNiV ’ ’ ‘*‘*‘’*' > 

‘*'*'*'f lu niKN 

I Ml 

I'M I It N lu u AL < Iluii A 

I’lM NMun auj » a, U ) 

I I t X t'Nl f /tKO» 01 J , MAlutk * 

N /naInx / Ndimi DlVSR, ObTI k« 

/IASI.,/ M.r.‘;,"i, "A.Kix.JS";,’''' 

/I MUX/'" Tmu iwavhsoj 

“-''’t'N /H ss,/ ; mm. 1 . iMiiIrioV 

/RAV.Nn/ M, •*• 

" ^ O 

>' < INO .m . ,1 , 

ULI |(J VVV9 

‘ m.7hcunuv rihuvu. 


■ *’ I '• 1 t »N(i|ni 
' n aii.m Am - 


- ALIM ) A, 


IS 


0IS.I9 

oisao 
OISaI 

OISA 4 
01 ‘»A A 
OJ bAf. 
01*'A#. 
oi;.A/ 
UlSAtl 
01‘*A V 
01 , '.;>(, 
OlASl 
OISA,* 

0 1 r> S 4 
Ol^AA 
OIS.S^ 

0 1 1> A A 

‘*l-AS7 

OlAAM 

0 I AAV 

01A6O 

OlAM 

OlAA/ 

OIAa.L 

01 AO A 
OlAoA 
OlAoo 

0 1 Ao J 

01 AOK 
OlAOV 
OlA/t, 

0 1 A 7 1 
OIA7.* 

0 1 A / 4 

OIA/4 
OlA/A 
01 A /O 
01 A /7 
OIA/h 

0 1 A /v 

01 AHO 
OlStt J 
01 A A.' 

01 ah t 

Ol'.tt., 

01 At. A 
UlAtlO 
01A|. / 

O 1 A A A 


III 


k 


o r> n r> r> r> 


I 
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